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INTRODUCTION 


The purpose of this study is to establish techniques for obtaining the 
minimum weight and power (or occasionally size) for a mass spectrometer 
system which meets a given performance requirement. Often this optimization 
procedure is necessary in the opposite direction where weight and pox^er are 
limited. The technique, however, usually proceeds equally x^ell in either 
direction, once the independent and dependent variables are identified. 

The type of mass spectrometer discussed in this study is the Quadrupole 
Hass Filter ^->2,3,4, and 5 w hich has been selected in several experiments 
for analysis of the earth* s upper atmosphere and potentially for planetary 
atmospheres. For the above purposes, it is typically operated in one of 
three general modes, as folloX'/s: 

a. Open Ion Source . Ambient particles are ionized xtfith limited 
energy (temperature) equilibration to the ion source; ionizing 
chamber gas density limited to ambient density as modified by 
partial accommodation to ion source temperature and surface 
conditions . 

b. Closed Ion Source, Close-Coupled to Ambient . Ambient particles 
are thermalized to ion source temperature; gas density established 

• by transpiration flow to ambient atmosphere or to a thermalizing 
chamber. 

c. Closed Ion Source, Excess Sample . Particles are thermalized to ion 
source temperature; gas density is established by flow through an 
inlet system from high pressure and is limited only by the leak 
rate as determined for the desired system operating conditions. 

For the purpose of this report a generalized configuration of a Mode 2 
ion source and the quadrupole is shown in Figure 1. 

Each mode of operation requires a slightly different treatment and will 
arrive at a different optimum instrument. For example, analysis of energetic 
particles with an open ion source will generally require a larger system in 
order to accept the diverse paths of the particles. Conversely, high ambient 
pressure generally permits use of a smaller instrument when the ion source 
can operate at its linear limit of gas density. Other modes of operation 
can arise xtfith a similar requirement for modification of the constraints. 
Hox<rever, this discussion is limited to the three listed. 



METHOD OF OPTIMIZATION 


Before an optimum instrument can be derived, a definition of the problem 
is typically given. For this study, since a range of applications was con- 
templated, attention was directed at maximizing the specific sensitivity of 
the mass spectrometer. Additional elements of the optimization procedure are 
required in order to arrive at a complete instrument. 

Expressing this analytically, the performance to be achieved can be given as 
total sensitivity: 


<s T ) 

i 



(e /n ) 
o o 


f 


= 1 


volts output/unit density, ambient; for q species, readable in the density 
ranges 


(n Q ) (min) < (n Q ) < (n Q ) (max), 
i i i 


under a given set of operating conditions and in a given environment. The 
solution to the problem requires a consideration of interfering species at 
various density levels, both directly as ion fragment current and indirectly 
as "noise" (e.g., neutrals and secondary effects); the system linearity at 
high ambient densities; the system noise at low ambient densities as well as 
the implied design and operational parameters. 


Assuming for each species a relationship between ambient and ion source 
gas density. 


n = H, n 
s 1 c 


( 1 ) 


where Hp is generally a time varying function relating both to the sampling 
mode and vehicle orientation. 

The current formed. Ip, by electron beam bombardment of the gas is 
substantially linear with density, 

it « K n . (2) 

F as 

+ 

The ion current collected I G can be expressed as, 

4 ■ s s 9 a 4 <3) 
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■where 0g = ion source transmission efficiency and 0^ = mass analyzer trans- 
mission efficiency, each of which is in general a function of the mass 
spectrometer operating conditions. 


For the purposes here, the output voltage, e 0 , can be assumed pro- 
portional to collected current, or 


e 

o 




c 


(4) 


combining (1) to (4), 


S„ = e /n = K K, 0 C 0. H . 
T o o a d S A 1 


(5) 


A knowledge of the variables under a given set of conditions, then 

allows determination of the unknown gas densities (n Q ) by observation of 

e . . i 

oi 


This study is then concerned with maximizing the product 


K 

a 


= S, 


the specific mass spectrometer sensitivy, in amperes ion current/unit density 
in the ion source, i.e., 


S = I /n = K 0„ 0 , 
c s a S A 


( 6 ) 


In this approach, all of the performance parameters are analytically 
related and summarized in the one parameter S which then becomes a function 
of size, weight and power. The optimum instrument is then inferred by the 
least value of weight (as a summary parameter) that satisfies the sensitivity 
and signal- to-noise ratio required in the problem. 

While is an important variable in the complete problem, it is only 
discussed to the extent that it impacts the mass spectrometer input con- 
dition in the mode of operation as described in the Introduction. 

The method of optimizing S takes one of two approaches, depending 
upon whether the quadrupole analyzer is operated in the "flat-topped peak" 
or "pointed peak" condition, (these terms are more fully defined on the 
definition pages of Appendix 4) . 

The first case where the instrument is designed and tuned to produce 
flat-topped peaks requires that ©a. - 1*0? i.e., 100% of the ion current 
received by the analyzer is transmitted and collected. The great importance 
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of this characteristic is that minor variations of the analyzer operating 
conditions or the particle initial conditions do not influence the accuracy 
of the measurement. This, along with an ion source 'that is to- first order 
independent of particle initial conditions (described later), is an essential 
factor in substantiating the measurements in certain space experiments. 

It is found that optimization can be approximated when 0^ = 1.0. Ana- 
lyzer entrance conditions of angle, position and energy of the particle will 
be solved in terms of desired resolution and peak shape. Secondly, this 
equation is substituted in an expression for 0s • 9s f s then differentiated 
for a maximum with respect to the remaining variables. 

In the pointed-peak mode 0^ (<1.0) becomes intimately a function of the 
electrical operating conditions of the quadrupole. Since there is as yet no 
closed, simple expression for transmission of the quadrupole as a function 
of the various conditions, optimization is a lengthly iterative procedure. 

In practice an open-ended solution of 0^ for a given set of conditions is 
usually possible. 

The quadrupole mass analyzer is first described in terms of its operat- 
ing' parameters in the next section. Subsequent sections, respectively, 
analyze the ion source parameters and combine the two equations. 


4 



ANALYZER CHARACTERISTICS 


By continuing the use of the optimization techniques of Reference 3 it 
is found that the quadrupole functions in a manner that appears similar in 
form to a double-focusing magnetic mass spectrometer with first order 
coefficients of angle and entry position equal to zero. It is thus an 
instrument of higher order. The equations can then be solved simultaneously 
with the ion source equations to optimize performance and establish 
sensitivity. 

In general terms, all mass spectrometers operate within limits imposed 
by the solid angle of ion source emission, , the area, A e , of ion emission 
and the upper and lower limits of an energy or velocity band, + 3 where 

3 ® — and v = velocity, 
o 

The quadrupole is no exception. It is the object of the optimization process 
to maximize this multi-dimensional "volume" to obtain the greatest current 
output. However, the quadrupole is different than most other mass spectrom- 
eters in the nature of imposition of penalties of weight and power in addi- 
tion to it being unique in providing versatility. 

Assume 'that the operation of the quadrupole analyzer is understood 
from prior work (c.f,, References 1-5). Further, the analysis assumes that 
the quadrupole functions in a theoretical manner with respect to ion en- 
trance conditions, i.e., the position, angle and energy of the ion as 
injected from the ion source arrive in a pure quadrupole field without end 
effects. Both quadrupole bxasing3 with nozzle injection and segmented rods^ 
attempt to nullify or improve upon the effects of quadrupole rod termina- 
tion; however, the full effects are not known. 


Scanning and Resolution 

Following common usage, the coordinate system is oriented with the z 
axis along the axis of the rod structure, y is through the negative pair 
of rods and x through the positive pair (see Figure 1). Then from the 
quadrupole operating equations. 
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8 e V 


dc 
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Cl) 


and 


q = 


4 e V 

ac 



2 

to 


(7) 
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where 


a, q 
e 



V 

ac 

m 

r 

o 

to 

a technique of 


= dc and ac stability parameters, 

= electronic charge (couloumbs), 

= applied dc rod voltage (volts), 

= applied ac rod voltage (volt-peak), 

= ion mass (kg.), 

= characteristic rod spacing from the axis (meters), 

= angular frequency of applied ac voltage (radians/ s ec ) , 
sweeping or scanning the spectrum can be described. 


With reference to the stability diagram in Figure 2- (a and b), which 
show a very select portion of the overall stability characteristics, the 
region of interest is around a = 0.23699 and q = 0.706. Charged particles 
are stable when 3y > o and o < B x < 1, consequently within the triangular 
region. Scan lines 1, 2, 3, 4, and 5 depict lines of constant ratio: 

a/q = 2 V dc /Y ac (8 


which pass through the stability diagram when a given mass-to-charge ratio 
m/e is in the quadrupole field and as V ac or a> is varied. Higher resolu- 
tion is achieved by increasing the scan line slope, a/q. Equation (7) can 
be used to establish parametric values which produce stability. 

In order to define the normal scanning performance, it is desirable to 
relate the peak shape to the stability diagram. Figure 3 depicts the ion 
current transmitted through the quadrupole as a function of q, (and as im- 
plied by constant a/q) as the peak is scanned. 

As a and q increase along scan line 3, for instance, ions of say 
(m/e) 1 are nominally stable in both x and y direction between A and D. 
However, a small, increasing current is transmitted and observed upon 
approaching A, which is due to ions of (m/e), unstable in the y direction 
which do not reach the rods for discharge within the limited length of the 
rod system. The longer the rods, the greater will be the fraction of 
unstable particles that reach amplitude r 0 and are eliminated. Brubaker 2 
shows typical paths while the Quadrupole Ion Entrance Mask study of this 
final report shows the envelope of motion of unstable particles. 

At the stability boundary at A, ions of (m/e), become y theoretically 
stable; however, they can also have a stable amplitude that is larger than 


6 



r Q depending on their initial angle and position off-axis. As the scan line 
is followed from A to B, ions with successively larger y angles and 
offset have stable amplitudes less than r 0 and are transmitted. If the 
angle and offset are bounded at small value, they all have maximum ampli- 
tudes less than r Q and 100% of the current is transmitted between B and C. 
The effects are approximately the same in reverse order upon crossing the 
x stability boundary from C to D and beyond. The correlation with peak 
shape is apparent from Figure 3. 

Figure 3 also shows the peak width characteristics defined as the 
following: 

Aq^, = peak top width measured in q 

Aq = theoretical peak width between boundaries measured in q. 

Aq B = peak base width to a defined fraction of peak height, 
measured in q (see definitions, Appendix 4) . 

Aq = increments of q taken by peak sides 
kx,y 

Aq = increments of q taken by tail. 

x,y 

These elements of peak shape will be used below. 

Considering, then, the technique of scanning, it becomes apparent that 
sweeping voltage has certain advantages. While the spectrum could be 
scanned by holding V ac constant and varying to especially if it is desired 
to continuously sweep a wide mass range, simple low power and precise 
methods of doing this have not yet been demonstrated. One advantage of 
the voltage scan is the simple time relationship that is easily produced. 

From (7) if a and r are constant, 

o 

m = K V (9) 

ac 

thus, if we let V = kt where t = time, 

9 ac 

m = Kkt (10) 

and the peaks are equally spaced in time. Voltage scanning is assumed in 
the calculations of this report. 
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It is also possible to maintain constant resolution during the scan. 
Resolution is defined in Appendix 4, but can be more readily visualized 
from Figure 4a. Two adjacent, equal sized peaks are considered resolved if 
their contribution to the valley between them is equal or less than a speci- 
fied fraction of the peak height. A stated resolution, m/Am is reached when 
the contributions are equal to the specified value. 


Taking this limiting case, two mass peaks are resolved when from (7) , 


3V 


ac 


V 


ac 


9V 


ac 


V 


m = cons tant 


ac 


q = constant 


thus, = 0 r, f or small increments, we assume 

q 


m 


Aq s 

— — = A m/m, where Aq g is the peak separation 


( 12 ) 


For the example in Figure 4a, 


Aq = Aq = q A m/m 
s a 


(13) 


In the more general case, the peak separation expressed by qAm/m can he 
greater or less than Aq B for adjustment of peak overlap, thus we let 

k Aq = Aq = q A m/m (14) 

Ifl. S 


For kjjj = 1, the peaks are resolved at the base by the conventional de- 
finition. For = 1/2, as shown in Figure 4b, the tail of adjacent peaks 
reach the measured peak center and the crosstalk becomes predominant. This 
we define by the tail fraction 0 at the edge of the peak base. Other degrees 

of overlap can be achieved by 1/2 < k <1. 

m 

In order to establish the conditions where each peak is just resolved 
from its neighbor (unit resolution) let 

A m/m = m /m or Aq,, ~ 0.706 m /k m (15) 

ii 4 R it m 


where m = unit mass, 
u 


From (8) 


a V 


V dc“ 


ac 


2q 


(16) 
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0.706’ 


and letting a = a , - Aa 
pk 

where = 0.23699 = value of a at the peak of the stability diagram and 


Aa^ yQg = increment of a between a ^ and the intercept of the scan 


line with q = 0.706, 


then. 


V 


V 


ac 


dc 2q 


(a P k - ia 0.706>' 


(17) 


However, Aa^ = 0.2524 Aq from the geometric relationship of the 
stability diagram and Aq^ — A^ + ^qt’ ( w ^ere 4q^ ~ Aq^ 
definition from figure 3) . 


Experimentally it is found that Aq t remains constant if mass and V ac 
are constant as resolution is adjusted by increasing or decreasing V^c* The 
peaks are seen to change in width in the central 100% transmission region 
without a marked change in shape or width of the tails. However, when 
scanning V ac and M, Aq t will change. From the Mask Study, Table I, Aq t 
appears to decrease with A m/m as m increases; i.e., approximately the 
same number of cycles are required for an ion to reach the rod at a given 
percentage of Aq outside the boundary whether the scan line slope is set 
for resolution 100 or 200. 


Thus, we assume 

(Aq t ) = k tl Aq, (18) 

where (Aq t ) is a specific tail width giving the desired quality of resolution. 

Also, for a given axial energy, the heavier ion will remain in the rod 
system for a greater number of cycles and thus have a higher probability of 
exclusion. Consequently, Aq^/Aq should even decrease at higher mass. 

Appendix 2 deals with the tail transmission probability and results in a 
relatively complex expression. For the purpose of describing the nature 
of unit resolution scanning we believe it is conservative albeit approxi- 
mate, to use (18) . 
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Thus, from (14), (17), and (18); Aq B = Aq (1 + kqq) 


and 


V dc = 


V ac %k 
2q 


0.2524 


V q m 
ac ^ u 


2q k (1 + k - ) 
^ m tl 


Substituting from the quadrupole equation 


4 e V 


q - 


ac 


2 2 * 
m r to 
o 


for V ac in the second term, and from fixed values for a ,/2q 


term: 


pk 


(19) 


( 20 ) 

in the first 


V, = 0.16785 V 
dc ac 


0.1262 


2 2 

q m r to 
_ u o 

4 e k (1 + k ) 
m ti. 


( 21 ) 


The second term thus expresses a constant offset in that is estab- 
lished relative to V ac at unity mass. 

Summarizing this section, it is shown that a continuous scan of the 
spectrum can be adjusted for unit resolution throughout the spectrum and pro 
vide that the masses occur at a uniform rate in time by letting V<j c = kq 
Vac - k 2 » and V ac = kt. Tor obtaining maximum information from such a scan, 
adjustment of k m allows for overlap of peak sides and tails such that a minx 
mum time is spent in measuring low utility information, k^q adjusts the 
base width allowance for the tails. The techniques of either stepping V ac 
from peak~to-peak or by compressed scanning are necessary to maximize the 
information content of a mass spectrometer. The quadrupole further exhibits 
the unique property among mass spectrometers of simple electrical adjustment 
to achieve both uniform pulse (peak) duration and repetition rate, as seen 
by the output measuring system, while maintaining a low ratio of lost time 
between peaks. 


Entrance Conditions for 0, =1.0 

A 

Tor 100% transmission through the quadrupole analyzer as considered 
here, it is necessary to limit the maximum entrance angles a x and ay; the 
initial off-axis position of the ion. 
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( 22 ) 


n 


X 


X 

e 

r 

o 


and n 

y 



o 


and the maximum injected energy from the ion source into the quadrupole. 

We establish these limits which govern sensitivity, through a tradeoff with 
resolution and peak shape. 

From Figure 3 we establish the base width 

Aq B = Aq + Aq tx + Aq fcy , (23) 

and the top width, 

Aq T = Aq - Aq x - Aq y , (24) 


where expressions for Aq^, Aq^, Aq^ and Aq^ are to be developed. 

Thus the "scan space" taken up by the sides and tails of the peak is 
given by. 


Aq B - Aq,p = Aq + Aq w + Aq^ - Aq + Aq^ + Aq 


l tx 


ty 


X 


or 


Aq T 

Aq B - AiT> “ Aq x + Aq y + iq tx + Aq ty. 

JJ 


(25) 


However, from (14), Aq B = q Am/kjn. Also, Aq T / Aq^ is the top-to-base 


m 


B 


ratio of the peak which we symbolize by T/B. Thus, 


Am 


, (1 - T/B) = Aq + Aq + Aq. + Aq. . 

km ' ' ^x ^y ^tx ^ty 

m 


(26) 


Consequently, we see that an improved peak shape (large T/B ratio) and 
higher resolution can be achieved by minimizing the Aq values on the right 
hand side of the equation. It is interesting to note that this is identical 
in form to a symmetrical magnetic sector instrument in which 

(1 - T/B) = 2 S ][ (27) 
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where Sp is the image or beam width and r the magnetic radius of curvature. 

Thus, the sum of Aq values of (26) is analogous to beam width. Compari- 
son of the dependencies of Aq upon entrance angle position and energy will 
disclose a further close resemblance of the quadrupole to the double-focusing 
magnetic mass spectrometer. 


For evaluation of Aq and Aq in terms of initial conditions we turn to 
reference 3 (Appendix l). y Their equation 14 restated here is 


|y max | _ 2.672 i 
.. Q 4= I 


0.505 r a 


JL 4- _hl 


B f, 

y i 


(28) 


where 


also 


max 


X 


V, 


a 


r 


x 


maximum amplitude of y motion 
initial (entrance) y position 
stability parameter shown in Figure 12 
(1. 0 + 0. 336 cos <j> ) 

0.336 sin <j> 

(V /V )l/2 s i n a . normalized transverse velocity 
I ac y 

phase of applied V at the time of ion entrance into the 
quadrupole field. 

voltage of ion injection into the quadrupole field 
initial angle of entrance in the y-z plane 


(V T /V ) 1/2 sin a 
I ac x 


This equation was derived directly from solutions to the Mathieu equa- 
tions describing the ion motion; it involved dropping higher order terms. We 
will use this equation for form and check coefficient values against later 
work. 


The largest value of attainable by given values of initial position 
y and normalized velocity y can be found by first letting |y max l = r Q - and 
solving for 3 , thus y 


y 


2.672 


,0.595 



(29) 
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where 


T) = y /r . 
y o o 

From Reference 3, the phase of fq and f2 giving the largest value to 
(29) occur at 


<i> 


3?r 


with little variation for a wide range of relative values of q and y . Thus, 
for this value of 4, Y y 


3 ~ 2.07 a + 1.09 n . 

y y y 


(30) 


We now let the maximum values (n ) and fy ’i determine (8 ) 

y max w y''max v y max. 


or to generalize, 

(6 ) = F- (y ) + F 0 (p ) 

y max 1 v 'y'max l y max 


(31) 


From the stability diagram we also note that due to the slope of the 
scan line in passing through lines of equal 8, 


Aq = 2.55 B 

y y 


(32a) 


and 


Aq = 1.30 (1 - 8 y 
x 


(32b) 


From computer evaluation^ of stable trajectories in the marginal trans- 
mission region (Aq x space) indicates that a function having similar propor- 
tionality to (31) exists for (1 - 8 X ) . Thus we assume 


(1 - 8 ) = F„ (y ) + F, (q ) 

^x y max 3 x / max 4^ x / max 


(33) 


Combining (26), (31), (32), and (33) 


q £2- (1 - T/B) = 2.55 (B ) 2 + 1.30 (1 - 8 ) 2 + Aqt + Aqt 

k m y max x max M x H y 

m J J 

(34) 

Letting the brackets ( ) denote maximum values for n and y , 

1 a - T/B) = 2.55 [Oq ( Y ) + F 2 (n )] 2 + 1.30 [f 3 (y j + f 4 („ x ] 2 
m y y 

+ Aqt x + Aqt y (35) 
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We would now like to simplify equation (35) and solve for the reduced 
coefficients in terms of some recent computer results 6 . This will allow us 
to update the coefficients to more reasonable results instead of using the 
theoretical values of (30) . We assume a relationship will exist for maxi- 
mum. transmission whereby 

(Y x ) = (y ) and (r^) = (n y ) and where k y 


and k. will have an optimum. Expanding and collecting terms we find that 

n 


q a - I/E) = G u (Y y ) 2 + 2 G 12 <Y )<1 ) + S 22 (Ey ) 2 + Aq tx + iq ty 
m 

(36) 


where 


G n = 2.55 


9 2 2 
F +1.30 F^ k y 


G 12 = 2.55 F 1 F 2 + 1.30 F 3 F^ k y k^ 
G 22 - 2.55 F 2 2 + 1.30 F 4 2 k 2 


In order to arrive at values for G;q, G l2 , and G22 which depend upon the 
density distribution of Y and p within their maximum values (among other 
factors) it is necessary at this point to rely upon the computer summation 
of a matrix of transmission probabilities in the stable region. For this 
comparison we have one value. In excess of 10 trajectories were run for 
a constant density distribution to (n x ) = (by) = 0.08, (Y x ) = ^y^ = 0*064 
and for each 0.1° of phase. Odd parity in the signs of n and Y give one 
transmission efficiency and even parity gives another efficiency. When 
averaged these values give the net efficiency. Sixteen values each of p and 
Y were run at the central value of each An = 0.005 and AY = 0.004. The 
quadrupole rods were assumed to be of infinite length so that no assumptions 
of transmission efficiency near the stability boundary or as a function of 
mass were necessary. Thus, Aq^ = Aq-j-y = o. 


The resolution intercept of 21 for 100% transmission in Figure 5 re- 
presents the case of T/B = 0 while resolution now pertains to separation at 
the theoretical base of the peak (k m = 1). Thus from (36), 



CY y V 


+ 2 e i2 <y <y + g 22 



(37) 
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We now note that this is the same result as if we had assumed from (35) 


that 


2.55 [F x ( Yy ) + F 2 (y] 2 = 1.30 [F 3 ( Yy ) + F^ (n y )] 2 (38) 

In fact we will find later, in Appendix 3, that this is a condition 
of maximum transmission for constant coefficients. It is logical that this 
would be so since the same q-space is taken up by the conditionally y-stable 
and x-stable sides of the peak. 

From (37) we define a perfect square 

G u <y 2 + 2 G i2 <y <y + g 22 <y 2 = (yy + ® 2 <y i 2 

and from (35) and (38) , 


= 2 x 2.55 [F 1 ( Yy ) + F 2 (n y )] 2 , (40) 

when T/B - Aq^ = Aq^_ = o and k = 1 

ty tx m 

It is further shown in Appendix 3 that the transmission will be maxi- 
mized when 

(y <Y y > = (6 2 )(n y ) (41) 

Consequently, from the computed intercept of m/Am = 21 at Y = 0.064 
and n = 0.08, 

(o.o64) y , 

and G = 1.43 = ( Gy 1/2 =2.26 F 1> 

1/2 

G 2 = 1.14 = (G 22 ) x/Z = 2.26 F 2 , 

G^ — 2.05, G-^ 2 = 1*64, G 2 2 = 1.315, 

F-^ = 0.634 and F 2 = 0.508. 

We could also have proceeded from (35) and (38) in the x-plane with 
2 x 1.30 [F 3 (y 4 ) + F 4 (n x )] 2 


0.706/21 


= ( G i 


(0.064) + G. (0.08) 


2 -( 


2 G, 


15 



to obtain 


F„ = 0.889 and F, = 0.709. 
3 4 


The above argument of coefficient values rests upon the choices of 
Y = 0.064 and r\ = 0.08 and the resulting resolution of 21 in the computer 
This is probably not the optimum since the choices were based upon 


run. 


preliminary results that emphasized the y-plane. However, to optimize the 
transmission data versus phase angle and n and Y for other Y/n ratios and 


W 


rj /p ratios would require an extensive effort. 

x y 


The choice, now, for an optimum could just as x-rell have assumed Yx and 
n x were (2 . 55/1.30) 1/ 2 larger than Yy and riy since this is the extra capacity 
in the x-plane. As a result F 3 and F 4 would then be of the same order as Fq 
and Fy» Since we have no logical basis for making this assumption we use 
the computer results with F 3 and F 4 values above. 


Returning to the comparison of coefficients from (40) ff, we find that 


(3 = 0.634 (Y ) + 0.508 (n ) 

7 y y 


(42) 


which yields substantially smaller coefficients than the theoretical values 
from Equation (30) where 


6 = 2.07 (Y ) + 1.09 (n )• 

y y y 

Smaller coefficients will lead to increased transmission since Y and n 
will be larger before limiting resolution is reached. It is important then 
to examine this difference. There appears to be two explanations, as follows 


(1) The G coefficients are in effect hybrid x and y plane coefficients, 
since the q space equations were made equal for the two planes. 

The scan line angle provided an advantage of 1.3/2.55 for the x 
plane, which when optimized reduced the effective y coefficients. 

(2) The computer provided a distribution of density versus Y and q. 
g obtained from (30) is the largest value reached by any ion; 
whereas in reality, we are dealing with a number distribution 
associated with the function F^ Yy + ^2 ^y which has the distri- 
bution shown below if ions are uniformly distributed in Y and n. 

Consequently, an apparently improved resolution is obtained at the 100% 
transmission intercept of Figure 5 because of the ion density at the limiting 
value . 
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CONST. 

CONST. 


Conseguently, we retain the more realistic values of (41) to describe 
the optimized (equalized x and y) 100% transmission edge. 

We return to (35) where a satisfactory solution has been obtained for 
the stable portion of the equation and now draw upon the results of Appen- 
dix 2 to describe the width of the peak tail. From (A2-21) approximate 
expressions for the q-space occupied by the tail to a given transmission 
fraction 0 are 



A V = 

0.0584 

n 

In 

0.0667 

0 

y 

0.00250, 


* 

0.0297 

n 

in 

0.0253 

0 

X 

0.00128, 

or if 

0=0, 
x y* 






Aq. = Aq. + Aq^ 
t H tx H ty 


0.088 ^ n 0.048 


n 


0 


- 0.00378: 


where n is the number of a-c cycles experienced by the ion and 0 X = fractional 
transmission at Aq^ x and Gy = fractional transmission at Aq^-y. 

f A 


Letting n - 




and combining with (42) yields 
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0 (1 - */»> - 2-55 [Fj (r y ) + F 2 (n y )] 2 + 1.30 [F 3 <Y x ) + F 4 (r, y )] 2 

0.88 (v ) 

+ in ^2. -.00378, (43) 

where f = frequency, Hz 

l = rod length, meters 

(v ) = maximum axial (z) ion velocity for the mass observed, 
meters/ sec. 

We now have an empirical equation defining resolution and peak shape as 
a function of initial position, angle and axial velocity. The top width is 
defined as the q-space taken up by the 100% transmission region of the peak 
and the base width as the q-space between equal x and y transmission frac- 
tions 0. The equation appears to be useful when the tail width, Aqt 
0.000706 or 10% of the theoretical peak width at a resolution of 100. 

This mass analyzer function must now be matched to the ion source 
function, which can also be described in terms of the maximum (emitted) ion 
position angle and axial velocity. 
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ION SOURCE CHARACTERISTICS 


It is desirable to update the ion source function in terms of particle 
density and the mode of use. 

The ion current formed in an ion source can be expressed as, 

* 4 * — 

I f = Q m 1 n s a e ^PS (44) 

where 

2 

= ionization cross section for ion of mass m (meter ) 

I = ionizing electron current (amps) 

3 

n = molecular density in the source (molecules /meter ) 

S 


= useful ionization path length (meter) 

If the electron beam has a useful width, w, for passage of the electron 
current I“> then the ionization density formed is. 


, Q I n & QIn 
j+ _ _m s_e _ s 

f wA w 

e 


which can be treated as an ion "cathode" with density proportional to elec- 
tron current per unit width (l"/w) . 

If magnification of this ion object takes place in the source indepen- 
dently in x and y, to form an image of the ion cathode at the final aperture 
as shown in Figure 14, then the illumination is, 

e w M M 
x y 

where and My are the respective magnifications. The dimensions of inter- 
vening apertures are assumed large enough so as not to remove current from 
the theoretical image. 
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The emitted current from the source is then. 


. > Q I n A 

, m s e 

I *= J A = T 7 — r: amps 

s e e w M M r 

X y 


where A e is the exit aperture area (meter ) 


( 46 ) 


When the pressure is limited only by instrument performance we use 
an. approach based upon space charge effects in the ionizing region. This 
case (Mode 3, excess sample) is considered first. 

7 

Brubaker has described the space charge conditions occurring in an 
ion source of infinite extent. His equations can be used to interpret the 
limits of pressure nonlinearity and therefore accuracy of the system. I - 
and n are then established by. 


and 


where 


I /w = 


i 

v R 
w 


amps 


n = X* = X 


4 d Q 


V m 
_r e_ 

a V mi 
e 


1/2 


(meter) 


(47) 

(48) 


i =» parameter for normalizing electronic charge density 
(amps/volt) 

V = Repeller voltage (volts) with respect to accelerator 
voltage 

3 

n = molecular density (molecules /meter ) in the ion source 
s 

X = fraction of critical density 

3 

n = critical density (molecules/meter ) 
c 

d = spacing of repeller to accelerator (meter) 

Q » cross section for ionization of ion mass m (meter^) 

m = mass of the electron (kg) 
e 
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a = fraction of distance of electron beam from the" 
accelerator to repeller 

= electron voltage (volts) 


m «= mass of ion (kg) 
We redefine 

M 

max 


Q 2 Q M"^ 2 m 2j ^ 2 (with M in AMU) 

m M = 1 111 u 

to sum all the charge present as intended, where m^ = unit mass (kg). 


Thus 


I' = 


Q A 
m e 

M M 
X y 


3 V, 


1/2 


w ^v V R? 


4 d 


IX ^ /2 



(49) 


The equation now shows the strong dependency on Vr reflecting that the 
rapid removal of ion charge provides a higher level of current for the 
source space charge conditions . For a nonlinearity (inaccuracy) the order 
of 5%, i v / w - 4 x 10“ ^ amps/volt meter and x - 0.6 from experimental data. 

In terms of pressure, for V R = 20 v and with a = 0. 5 , V e = 75 v, d = .0025 
meter, the critical pressure is 3.6 x 10“ 4 torr. Thus when x = 0.6, the 
operating pressure can be 2.2 x 10“4 torr while maintaining linear operation 
within approximately 5%. 


Gathering terms, 

K A V ^ 2 
+ _ 1 e R 

s M M d 
x y 


amps . 


(50) 


where 


K l = 


3 1 x Q 
v A m 


4 w/^Q M 


1/2 


m 


1/2 


a V m 
e u 


amps 


(volt) 3/2 meter 


Assuming that the ion source will be used with a quadrupole mass ana- 
lyzer it is desirable to express the strong dependence of (50) upon Vr in 
terms of the axial velocity through the quadrupole, v z . It is noted that 
there is a tradeoff between the added sensitivity to be gained by increasing 
repeller voltage and the loss in resolution from increased ion velocity 
spread and thus larger tails. In the limiting case, the quadrupole axis 
can be biased to the potential of the least energetic ion created in the 
electron beam without losing ions by reflection from the quadrupole 
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potentials. Thus the ion of maximum velocity (created at the most positive 
potential of the electron beam) will have an axial potential. 


V z “ d V R 

where t is the thickness of the electron beam. 

In practice it is often desirable to have a slightly larger v z than 
described by the real electron beam thickness. However, one may describe t 
as a pseudo- thickness which remains proportional to Vr. Thus, 

m CvJ 2 

t V R /d = jg volts , 


V R = 


dm (v ) 
z 


volts , 


and (50) becomes 


Tr . ,1/2 t s 3 

T+ K 1 A e d 

s MM 

x y 


(*) 


Q 

Further Abbes' Law can be shown to describe M and M as follows: 

x y 


sin (a x ) , 


sin (a ) . 

y 


where V Q is the initial transverse ion energy per charge to be accepted by 
the ion source lens system and (ot x ) and (ay) the maximum emitted angles of 
the ion source for the value V Q , 


Substituting in (51) and letting 

\ 3/2 rt I 


- ! 1/2 (*) 


3/2 3 4 

volt sec /meter , 


+ _ K 1 K 2 A e V I V z 
T s " V 

o 


sin (a ) sin (a ) amps 
x y 


22 



Preparing for combination with the mass analyzer we let A e = 4(x e ) (y e ) 
“ ^(nx)Cny) r 0 2 from (22), and Vx sin (a x ) sin (ay) = Vac (y x ) (y v ) from 
(28) . Thus : y 



s 


4 LL V r 2 

v c 0 (v z > amps 


(55) 


where again ( ) denotes the maximum value of the included parameter. 

Equation (55) now expresses the ion source emitted current as a function 
of variables common to the mass analyzer. It demonstrates that increased ion 
source current can be obtained by a larger and more powerful quadrupole or 
by increased ion velocity, angle or entrance area. 


If we now consider the ion source equation when an excess of sample is 
not available we find that ion space charge is no longer a limit on perform- 
ance. We retrace to Equations (47) and (48) and accept Equation (47) as valid 
since it describes the space charge depression due to the electron beam. 

The limiting ionizing current will be reached when space charge produces 
the same potential in the ionizing space as at the accelerator; i.e., zero 
field gradient in front of the electron beam. The potential then drops from 
Vr at the repeller to zero at the electron beam creating a gradient, 

E ~ V * 

Z " (1 - a) d 


in this space. From Poissons equation, 


d 



- - p/e 


o 



or 


in the infinite planar model assumed. 

= potential in the ionizing space (volts) , 

3 

p = spare charge density (coul/m ) 

£ q = permittivity of free space. 


The gradient change through the electron beam is then 
d E Z 

= P /e = J /v e 
o e e o 
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where 


J = electron, heam current density amps/meter , 
e 


and 


v - 
e 


2 e V 


m 


1/2 


electron velocity meters/sec. 


The useful electron current passes through an area of wt where t is the 
electron beam, thickness and w the useful width. Thus 


AE = 


I Az 


z w t v e 
e o 


(56) 


But Az = t. Consequently, A E z = I /w v g £ q or, for Modes 1 and 2, 


v e V_ 
e o R 


(I /w) l,2 " V e £ o A E z ~ (1 - a) d 


2 e V \ 1/2 e V 
el o R 


m 


(1 - a)d 


amps /meter 
(57) 


Equation (48) does not apply since there is assumed to be negligible 
ion space charge. In practice, of course, either positive charge would 
build up or slightly less electron current would be injected to create a 
slight gradient toward the accelerator. 


Substituting (57) into (46) gives 


X = 


Q e 

hn o R 


2 e V 


1/2 


s (1 - a)d y 


n A amps 

S 6 


(58) 


Proceeding along previous lines and substituting (51), (53) and the 
equations following (54) into (58) we have 

. 4 1L n V r 2 „ 

I*- 1 2 y s . ac -g- <v z ) (n x ) (n y ) (Y x ) (Y y ) amps (59) 


where 


and 


K. 


i Q e 
' o 


2 e V 


1/2 


1 (1 - a) 


m 


amp meter /volt 


| ft O 

K„ = -zrz volt sec /meter 
i zet 
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In general. Equation (59) is used whenever the ambient density limits the 
ion source density, n s (either Mode 1 or 2) , and Equation (55) when excess 
sample is available (Mode 3) and thus positive ion space charge limits the 
source pressure. Mode 1 then differs from Modes 2 and 3 in use of an 
imaging ion source as described in the above equations. Modes 2 and 3 
could, in general, use a thermally imaging ion source since the initial 
energy of the ions is controlled by the ion source temperature and is not 
excessive. We do not go through the exercise here of describing the 
thermal-focus source because it gives substantially the same end result and 
requires a more emperical approach. 
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OPTIMIZED INSTRUMENT FUNCTION 


The maximum exit conditions of the ion source from either (55), (Mode 3) 
or (59) , (Modes 1 and 2) now provide the maximum acceptable input conditions 
for the mass analyzer as given in (43) . Since we have predicated the devel- 
opment of the mass analyzer equation on 100% transmission, the current trans- 
mitted by the ion source becomes the collected current to the detector, i.e., 
I 4 * = I+. Thus, 

|"H ■* 


and 


4 LE V r 

C c = v £—2 ' ( 0 (V <V CO <Y V ) amps, 

3 o y y 


(60) 


and 


' 1,2 


4 k'k" n V r 2 

s_ac — o ( v (0(0 (0(0 am P s > 


V 


x y ‘x y 


(61) 


where 


K l = 


3 i X Q 
v m 


m 


1/2 


4wEQ m H' 


1/2 \ a V m 
1 e u 


amps 
. 3/2 

volt meter 


.1/2 ( m \ 3 ^ 2 volt 3 ^ 2 sec 3 

K z " IIF j ~h — 

' meter 


K ’ - £ ° 
1 " (1-a) 


2 e V 


1/2 


e 


m 


amp-meter 
volt * 


' m volt sec 

— 


"2 2et 3 

meter 


In Appendix 3, Equation (43) is combined with the above ion source 
equations by solving for (y x ) , and substituting in (60) and (61) and then 
differentiating for a maximum as a function of (n x ) > (by) , (ly) and (v z ) . 
In this process the optimized parameters are assumed to be independent of 
each other and that the coefficients are constant. 
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The resulting maximized Mode 3 (closed source) equation (A3-22) , is, 


0.957 Kl K 2 V ac CH) 3 r 2 

I (max) = 

C 3 


f iW« v 


(t .048 \ 

<M n —j 


!(v( 1 _ b) + - 00378 ) 5 ^ 


(62) 


where optimized values are. 


< v z >3 


( V3 


1 (f&)(U+ .00378) 

5 .088 In (^M) 

2 (U + .00378) 

2 P, | 5 (2.6) 


1/2 


(63) 


<V>3 


2 F, 


2 L 


2 (U + .00378) 
5 (5.1) 


1/2 


( V 3 


2 F 


3 L 


2 (U + .00378) 
5 (2.6) 


1/2 




2 F- 


2 (U + .00378) 
5 (5.1) 


! 1/2 


and 


U = 


q&m 


k m 
m 


l 1 ' b ) • 


For Modes 1 and 2 (open source) from (A3- 26) 


<maI °l,2 


0.152 K*k' n V (f£) 2 r 2 , A 
12 s ac o / qAm 

2 l k m 
m 


F F F F V 
1 2*3*4 




(“!)*• 


00378) amps 
(64) 


The optimized values of the parameters are somewhat different due to 
the power of (v ) ■ Thus , 


, , 1 (f£) (U + .00378) 

z1 ’ 2 = 2 .088 in ’ 


(65) 
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^ n aPl,2 


( Vl»2 


( ' Y x ) 1 s 2 


C Vl»2 


2 F, 


2 J 


2 F„ 


2 F„ 


(TJ + ,00378) 

2 ( 2 . 6 ) 


2 L 


(U + .00378) 
2 (5.1) 


(U + .00378) 

2 ( 2 . 6 ) 


1/2 
1/2 

1 1/2 

1/2 


( U + .0037 8) 
2 (5.1) 


Discussion of Parameters 


The important conclusions at maximum transmission are as follows: 

a. The q-space taken by the ri and y functions are shared equally in 
a given plane, i.e. , Fi(y y ) = F 2 (n y ) and F 3 (y x ) = F 4 (n x ) • This 
is so, regardless of future improvement in the optimum values 
found for the F functions . 


b. The q-space taken by the x and y functions are equally shared, 
i.e., 2.55 [Fi(y y ) + F 2 (n y )] 2 =1-3 [F 3 (y x ) + F 4 (n x ) ] 2 . This 
results also in the conclusion that 2,55 gy-s - 1.3 (1 - 3 x )g at 
the edges of the 100% transmission band. Thus, (1 - 3 X ) E = 1.4 gy E 
where the subscript E now refers to the resulting practical edge 

of the 100% transmission band. 

It is noted that as a result of Conclusions a and b, the q-space 
by each of the four functions, Fi(y y ) , F 2 (n y ) , F 3 (y x ) and F 4 (r) x ) 
are all equal. 

c. The leverage in v z exerted by the "tail" function which is now 
included in the optimization equation is quite large. This effect 
is a result of being able to increase ion current for the same 
ionizing space charge when an increased repeller-accelerator 
voltage, and thus axial ion velocity, is allowed. There are 
several practical limitations to extending this advantage very far. 
However, one important influence is the election beam thickness 
which also governs the AV to be accepted by the analyzer. A well 
focused electron beam, such as from a line filament, can have the 
same importance as a considerable increase in rod length. 
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Conclusion c is only one of a series of parametric effects on the 
equation which can lead to practical improvements in transmission. We now 
reconstruct the equation in order to understand the physical effect of the 
various parameters. 

We start with Equation (60) and recall that the area of the exit aper- 
ture of the ion source is Ag = 4 r§ n x q y . Calculating from (63) for (n x )3 
and (hy )3 we find from our choices of F 4 = .709 and 1?2 = *508 that n x = n y 
or that we have a square exit aperture. (It is also recalled that an 
alternative assumption would have made F 2 «F 4 an ^ therefore n x r| y *) If, for 
example, we let Am/k^ = 1/50, T/B = 0.5 and with q = 0.706, U = qAm/k m m 
(1 - T/B) = .00706 then (nx)3 ~ (n y ) 3 = 2.88 x 10 - 2. Thus for an r 0 of 
2.5 mm, the exit aperture is a square, 0.144 mm on a side. 

Further, we recall that 


V 

ac 


Y Y 
x'y 



sin 01 sin a 
x y 


1 


M M 
X 


Calculating again from (63) and the above values (Yx ) 3 " 
For Vj = V ac - 208 V and V Q = 0.208 V, then sin a x - sin 
(a = 1.3°) and 


(Y y ) 3 - 2.3 x 10 2 . 
ay = 2.3 x 10“^ 


2 . 

M 



h - 0.727. 
M. 


Thus a small magnification of 1.37 is required between the ion "cathode" 
(the ionizing region) and the source exit slit. If we demagnify the exit 
dimensions , we find that ions need' only be drawn from an area 

(A ” (0.105 mm ) 2 = 1.095 x 10~ 8 M 2 

0 3 M 2 


in the ionizing region. 

In fact, of course, the first ion source aperture is made larger to 
assure uniform focusing properties in the center of this lens. In addition 
the focusing model assumed an imaging ion source (see Appendix 1) . For the 
more commonly used thermal- imaging source a larger first aperture would be 
used to focus the center of the thermal image into the exit aperture. 

We have, so far, found that 
2 

4 V r 

— fr — ° (n ) (n ) (y )(y ) = a 

V x y x y o 

o J 

inhere Aq is defined as the area in the ionizing region from which ions are 
drawn to achieve the optimized transmission. But, by the same token, the 
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collected current is then the ion cathode density times the effective 
useful area or 1^ = Jp Aq. Thus from (60), we must have that K 4 K 2 (v z ) ^ 

- < 4 )° ' ■ ' ■ ' " 

equations. Likewise, 


vhich could also have been shown by returning to the original 


<4>1.2 “ K i K 2 ", <V 2 ' 


Consequently , 


(<) 3 = [#3 IA 0 ] 3 - 1*3*3 (v 2 ) 3 ] 


4 V r2 

-• f ' ' ° ( V C V < V V 


( 66 ) 


and 


^ X c^l,2 '•Vl 5 2^ A o^l ) 2 ^ K 1 K 2 n s ( - V z- > ^ 


^f^(n x )(n y )(y x )(Y y ) J 


(67) 


(66) and (67) can be further expanded by the substitution of their 
respective optimized elements from (63) and (65). However, since the 
physical meaning is clear in (66) and (67) we retain them and turn to an 
examination of their first terms. Looking at the definitions ff (61) we 
see that none of the coefficients in Kj can have much leverage or sensi- 
tivity. The i v x/w product sets the nonlinearity vs density to be accepted 
and cannot be increased greatly without reaching space-charge effected ion 
flow. Operation is also generally near maximum cross section, Q m , by choice 
of V e . a has only small leverage. 

In K 2 now, we find that the election beam thickness, t, can have a sub- 
stantial leverage on sensitivity, further, the ion mass appears. The same 
is true on K^. We therefore isolate t and separate M and m u in the constants. 

Since M appears in several places it is best to note their differences 
noxtf. T.QJjin appears in the denominator of Kp. This reflects the total 
ion space charge summed over the spectrum while M reflects the velocity 
in moving toward the accelerator electrode. If we had traces of various 
gases, to say 100 AMU, in air the bulk of the charge could still reside in 
the low mass region and would remain constant. This is the 

assumption that is carried through the remainder of the discussion. M also 
appears in Equation (51) for and in K 2 end K This pertains to the 
mass in view (being scanned) as it moves in the repeller-accelerator field. 
Later on we will substitute for V ac and find M describing the mass in view 
in the quadrupole. The latter two are thus the same. 

Kp appears to have leverage by making a approach 1. While having the 
electron beam close to the repeller increases its capacity to the sur- 
roundings and thereby reduces the space charge depression, it also leads to 
practical problems of bombarding the repeller. It is felt that in practice, 
only limited improvement can be obtained here. 
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The constants are now calculated assuming the following values for 
illustration: 


_3 

i /w « 4 x 10 amps/volt meter 

v - 0.6 


Approximately 5% + 
nonlinearity in I vs r] 


That an sample is used having 


^28 


-20 2 
1.8 x 10 meter 


-20 2 

Q- . ~ 0.9 x 10 meter 
14 


£ Q m M 1 ^ 2 » 1.8 x 10 20 (28) 1/2 + 0.9 x 10 20 (14) 1/2 

= 1 , A -20 . 2 . ,.„ T v 1/ 2 

« 12.9 x 10 meter (AMU) 


a 

— 

0.5 

V 

e 

= 

90 volts 

( 


\ ^‘ //2 -3 -1/2 

\ » 3.49 x 10 (volt) ' 

V a V e 

m 

u 

) 

d 

J=3 

-3 

2.5 x 10 meter 

m/e 

~ 

28 x 1.035 x 10 8 kg/coul. 

K 1 

~ 

-7 3/2 

8.8 x 10 amp/volt meter 

K 2 

~ 

1.87 x 10" 10 t" 3 ^ 2 M 3/2 volt 3/?2 sec 2 /meter 

4 

“ 

„ „ „ -24 2, . 

1.8 x 10 amp meter /volt 

K 2 

= 

5.19 x 10 3 t ^ M volt sec 2 /meter 3 

K 1 K 2 

= 

, , a -20 -3/2 „3/2 2 , , 5 

1.65 x 10 t M amp sec /meter 

hh 

= 

9.28 x 10 33 t ^ M amp sec 2 /meter 
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Letting f = 4.0 MHz, l = 0.1 meter, 0 = 10 ^ and U = 0.00706 as above we 
can calculate the velocities from (63) and (65) . Thus 

(v ) = 4.79 x 10 3 meter/sec (m/e 28 at 3.2 V) 

Z j 

( v ) - 0 = 3.99 x 10 3 meter/sec (m/ e 28 at 2.22 V) 
z i , 2 

Also if xte assume n s for Mode 1 and 2 to be 2.42 x 10^/m 3 (10 ^ torr at 
400°k) and t = 5 x 10“^ meter, we find for M = 28 that; 

+ 3-2 2 

(Jj.)^ = K 1 K 2 ^ v z^ 3 = 2.38 x 10 amp /meter , 

and 

CJ lh,2 " K i K 2 n s ( Vl.l = 2 -° x 10 ' 2 “ p/meter2 - 


It is interesting to note that the formed ion current densities are nearly 
the same, yet the imposed conditions were substantially different. We go 
on to evaluate CAq) x , 2 f° r t ^ ie satIie values as before and find 

■(A q ) i 2 = 1.713 x 10 8 meter 3 . 

+ 

Combining to find I c for the above conditions, 


(I^) 3 = C J y) 3 C A 0 ) 3 = 2.38 x 10~ 2 x 1.095 x 10 8 = 2.61 x 10“ 10 amps 


(I + ), _ = (it), n (A ) n = 2.0 x 10“ 2 x 1.713 x 10' 8 = 3.426 x 10 
c 1,2 F 1,2 o 1,2 amps 


( 68 ) 
-10 


(69) 


Again, it is surprising to find the values so close. Looking for an 
explanation, it is apparent that the A Q values will be nearly the same since 
their optimized components are about equal. However, the arbitrary choice 
of n s for the open source should have provided ample difference in Jp. 

If we break J* into 
F 



we find that the distinguishing difference between the too modes is that 
the open source (Modes 1 and 2) has a high electron beam density but a low 
particle density. The reverse is true for the closed source. 
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From (57), 


(I /w)^ 2 = 0*459 a/m = 459 ya/mm, (70) 

and by arbitrary choice above, 

IQ Q 

(n ) ^ 2.42 x 10 /meter . (71) 

S aL } M 

Whereas from (47) , (51) and (63) 

= 0.0667 a/m = 66.7 ya/mm (72 

3 

and from (48, (51) and (63) , 

(n ) » 0.6 n = 0.6 x 3.32 x 10 19 = 2.0 x 10 19 /m 3 (73) 

S o c 

So it appears that a trade off of extremes resulted from the choices 
made. If we examine the practical validity of these choices, the current 
density of (70) probably cannot be conducted^ through a small isolating, 
rectangular slit if that is a requirement. Without setting up a whole new 
parametric system it is assumed that approximately 200 ya/mm could be used 
in practice if the passage length did not exceed 1 mm. Even this large 
value does point up that considerable sensitivity is to be gained by higher 
electron currents than the 10 - 40 ya range typically used. Sensitivities 
are therefore reduced by 200/459. Thus the revised values are: 



(I /w)^ 2 = 0.200 a/m, 

+ -10 

(I c ) = 1.5 x 10 amps, 

4- -3 2 

(Jp)l 2 = 8.8 x 10 amp/meter , 

K ^2 = 4.03 x 10 33 t ^ M amp sec^/meter^, 


(7QR) 

(69R) 


and 


f 



0.435 Q e 
m o 

(1 " a) 



(61R) 


18 3 

The density in Modes 1 and 2 of 2.42 x 10 /m can probably be sus- 
tained without serious nonlinearity ; however, it is not practical to say 
for sure without an actual test. 
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Now for the density in Mode 3. If a closed ion source is operated at 
the density of (73) , corresponding to 8.3 x 10“^ torr at 400 °k, the gas flow 
from an ion source of 50 cc/sec speed would be 4.14 x 10”^ torr cc/sec. 

This then, would produce a pressure of 4 x 10 torr in a 10 liter/sec ion 
pump. 


Without again optimizing an undefined task, it would seem that this 
sample flow and the resulting constraint is unrealistic for a small space- 
borne instrument. The limiting density is therefore arbitrarily reduced by 
a factor of 4.0 to 

(n ) = 0.15 n = 5 x 10 18 /m 3 (73R) 

s J c 

by making y = 0.15 and thus 


K 1 

= 2 . 2 x 10 ^ amp 


K 1 K 2 

, ,, ,.-21 -3/2 m 3/2 2 / 5 

= 4.12 x 10 t M amp sec /meter 


(<)s 

- 6.52 x 10 amps 

(68R) 


and 

= 5.95 x 10 3 amp/meter 3 

In principle it is possible to increase i v if x is decreased in order 
to maintain the same space charge and nonlinearity. However this requires 
an experimental maneuver since the theory does not provide the means of 
making the translation from one to the other. In addition, other options 
of increasing electron current density from 66.7 pa/mm while maintaining 
low ion source pumping speed or of enclosing the filament are involved. 
Since this additional complexity is not warranted here, the revised values 
are used as attainable numbers for illustrative purposes. 

The net result of these alternatives is to make the closed and open 
ion sources look more nearly alike in their operating levels. The fact 
remains, of course, that the size of the tube or leak to the environment 
is greatly different resulting in different source time constants. Because 
of the similarity of end results, however, only the open source modes are 
plotted; remembering also that for this source a larger electron density 
than typical is required to make the ion current values attainable. 
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Parametric Performance Gurves 

Using the Mode 1 and 2 equations with, the revisions discussed above, 
we proceed to plot the collectable ion current as a function of the various 
parameters. 

Retaining the variable elements of (64) we recall that 


(max) l,2 a 


M V (f£) 2 r 2 . . w 

ac — ( SM (! _ t/B) +.00378 ) 4 . 


‘ v o (/» "T) 


2km 
\ m 


(74) 


/ 


In order to reduce the number of variables we let M = and express 
the quadrupole equation as 


2 2 

V = K M r f 
ac v m o 


(75) 


where 


K = 7.2 x 10 

v 


M = maximum mass to be scanned (AMU) 
m 


f = frequency (Hz) 
and substitute in (74) to obtain 


+ 

I c (max) x 2 a 


,4 n 2 4 „2 

f A r M 
o m 


t V 


(u^f) 


qAm 

k m 
m 


(1 - |)+ .00378 


)‘ 


(76) 


Since power is also an important function we let 


P = 


2 V to C . 

1 = K V 2 fZ 

Q p ac 


(77) 


where Cq is the quadrupole rod capacity (farads) which is assumed to scale 
linearly with Z by neglecting end and fixed capacities . Q is the circuit 
quality characteristic tjhere it is assumed that by clever design Q can be 
held constant. While neither of these assumptions is valid in practice, 
the added complexity for a more complete description is beyond the scope of 
this discussion, 
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It is further an important simplification that we let £ be propor- 
tional to r Q for most of the description. To do otherwise is to invite a 
lengthly chain reaction influencing the axial velocity and the repeller 
voltage where such is not necessary. The independent effect of £ will 
be shown later. Thus for this relationship we use 



Scaling Size at constant power, 9 as a parameter .- The simplest 
variation that can he found is by letting; 

P = const., fr = 4 x 10^ x .00254 = 1.016 x 10^ Hz, meter 

’ o 

-4 

- const.. M = 50 const, t = 5 x 10 m = const, 

* m * ? 

V = 0.208 V = const, Am/k m = 1/50 = const, 

o m 

kjj, = 1, T/B = 0.5 = const and V &c = 208 V. 

+ 2 

Scaling ion current with size we have I a r . Figure 6 shows the col- 
lected current at P = 10“^ torr (400°K)and sensitivity as a function of 
r 0 with 0 as a plrameter. Sensitivity as used throughout this section 
is considered to be the same as m/e 28 from N 2 even though mass is freely 
changed to describe the instrument performance. 

It should also be noted that V Q = 0.208 V is the voltage of a 6 kT 

particle at 400°K. This is slightly more energy than required to encom- 
pass 99% of the one-directional velocities. 

Scaling size, H, as a parameter .- Here we deal only with the maximum 
mass to be scanned, M^, since this is generally the condition of maximum 
constraint. We further assume that to reach the maximum mass Vac and V^ c 
are increased with a fixed offset in to produce unity resolution, i.e. 

Am/m = 1/Mjn (81) 
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In order to set values for kjjj and 0 , we assume that the peaks sides 
overlap such that the tail transmission, 6, is just reached at the top 
corner of the adjacent peak. Thus the valley between the peaks is essen- 
tially filled yet the full top width is available for value determination. 


When this condition is reached, i.e., having an appearance lying between 
Figure 4 a and 4b, then 


qAm 

m 





or since 


and 


“ kill fr ° m C14> ’ 

m 

K - k ( 1 + i) > 
qAm ( 2 ~ I) 

0.5 m (l + §) ' 


(82) 


(83) 


Substituting (81), (79) becomes 



c 


(max) 


a 




(84) 


We now let : 


and 


fr 

o 

= 

1.016 x 10 Hz = C 
m 

r 

0 

= 

variable 

M 

m 

= 

parameter = 25, 50, 100, 

t 

= 

5 x 10 _4 ■» C 

V 

o 

= 

0.208 V = C 

0 

= 

_5 

10 

T/B 

= 

0.5 

V ac 

a 

M . 

m 


200 AMU 
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2 2 

Thus power scales as V ac or 1^. 

At r G = .00254 m and 1% = 50, V ac = 208 V, P = P 0 and I~/w = 200 ya/mm, 
= 1.74 x 10 ^ amps 

Again sensitivities for both I /w = 200 ya/mm and 20 ya/mm are plotted 
vs r Q in Figure 7. 

It can be noted that the sensitivity does not fall off rapidly at 
M m = 100 and in fact recovers slightly in going to n 200. This is the 
effect of the constant .00378 (associated with U) in combination with 
increasing V ac . This does not seem normal or logical and probably indi- 
cates that the narrowed peak width is within 0 . lAq of the theoretical 
boundary at a resolution of 100 and that the data in this range should not 
be used. 

Additional study of the tail equation will have to be done before this 
can be clarified. The next curve shows this even more strongly. 

Sensitivity vs Mass and Resolution .- This set of conditions is the 
same as in the preceding paragraphs; however, the role of M m and r Q are 
reversed. Figure 8 shows N 2 sensitivity vs with r 0 as a parameter. 

V ac is scanned to cover the mass range to M m while holding unit resolution 
as defined above. 

r Q = 1.27, 2.54, 5,08 and 10.16 mm and I /w = 20 ya/mm. 

It is apparent in looking at Figure 8 that the equations in their 

present form should be used at M < 100. 

m — 

Constant Size, Constant Power Stepping .- If we now wish to cover a wide 
mass spectrum by stepping mass range on a given instrument, several ap- 
proaches are possible. In each case V ac is scanned within the range and 
only the locus of sensitivity values at the top of the ranges (V ac (max)) 
is shown. 

One approach yields constant power by letting V ac (max) = constant and 
adding external capacity to hold = constant. Thus with t, V 0 , r Q and 

0 constant . 

1+ (rax) a (l - |) + .00378^ (85) 

and P = V ^ fc = constant 

,P ac 


I (max) 
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Inductance L = 1/w C = kC. 


Using k - 1/2 (1 + T/B), Am/m = 1/M and T/B = 0.5, 

m in 


+ .00378 


with 


r c (max) a UM 


r = .00254, I /w = 20 pa/mm, t = 5 x 10 , V = 0.208 V, 

o ° 

6 = 10” 5 and f = 4.0 x 10 6 Hz at M = 50, the sensitivity 


is plotted vs 3^ in Figure 9. 

Another approach is to scale V ac upward as the maximum mass in the 
range is increased. This reduces somewhat the steep fall-off in sensitivity 
with mass. We let C = constant, L a N 2 >V ac a I^/N^and f a (LC) 1 ' 2 

a 1/^2 where ^ and Ny are the secondary and primary turns of the r.f. 

transformer. Thus P a (Ng/lS^) 2 (N^/N^ ) = ^/^ a V^; M 2 f 4 a V & 2 a P 2 

a (^/N^) 2 and from (80), ^ a V^/f 2 a ^ (l/^) 2 “ N 2^1* T1:ius 

N„ a IL 1 / 3 M 1 / 3 . If we hold N- constant and substitute 
2 1m 1 

f 4 a 1/(N 2 ) 4 a l/4 /3 . 


I + (max) 
c 


a M 2 ^ 3 


Q 

1.5 M 


+ .00378 


using the same values as for (86) . The sensitivity is also plotted in 

Figure 9, using S = 1.74 x 10“ 7 a/torr at = 50 and I /w = "20 ya/mm. 

t, V and 6 are constant as above, 
o 

It is clear that without increasing size or by applying considerable 
power, sensitivity drops steeply with the maximum mass to be scanned in a 
given range (or conversely increasing with lower mass) . Only by techniques 
such as scanning V ac directly with M for the full range can the sensitivity 
be maintained. 

One thing that should be remembered is whenever resolution, Am/m, 
is changed in ^ ^ ^1 — — ^ + .00378^ that this implies a change in n. 


+ .00378 


that this implies a change in n. 


y and from the appropriate transforms, y automatically changes if V ac is 
the independent variable to cause Am/m to vary, v is partially corrected 
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and also 


as the mass in view changes , since v z 


1 

2 


(£A) (U + .00378) 

non 4 .048 

.088 In — 


(2et V E 



1/2 


Thus some compensation is made as m and U change but not by the same degree. 
This residual v z ,plus q then must be changed by some electrical or mechani- 
cal alteration in the system in order to achieve the sensitivities shown. 

One possible means to automatically reach the resolutions called for is to 
decrease Vr as mass increases. 


Effect of Length and Electron Beam Thickness on Crosstalk .- If we look 
at the effect of length and electron beam thickness on crosstalk, the first 
step is to hold the sensitivity constant. By also letting f, r OJ Mgj, V Q , 
Am/m, k^ and T/B be constant and reconstructing Z ~ kr^, which contributed 
r^ to the equation, 

then y2 

— - = constant (88) 

t (fn 

m and 9 = 10 k = 0.279 and 

(89) 

l with t as a parameter. 


For values of t = 5 x 10 m, & = 0.1 
6 = 0.048 exp - (&/(kt) 1/2 
Figure 10 express, 8 as a function of 



It is apparent that the effort to reduce t is well worth while if the 
length is an important physical constraint. 

Open Ion Source, Mode 1 .- For the focusing of energetic particles as 
is implied by the use of an open ion source it is necessary to accept ener- 
gies eV c that are substantial. The scaling of a system to accommodate 
energies of the order of 10 ev can be readily done with the optimization 
equation expressed in its present form. The principle impact on design is 
not only that system will be larger but that an ion source having a different 
character of focusing is needed, i.e., a true imaging lens system as des- 
cribed in Appendix 1. The present equations describe this type of source. 

Several scaling factors are simple to show if we plot vs r Q with V D 
as a parameter. For instance, we can hold power constant (at constant) 

giving I a r^/Vo. Or we can let f = constant to yield I c a r£/V 0 whereby 
5 C 

P a r G f, greatly expanding the power consumption. A compromise might be 
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f 2 = k/r Q in which P a r 3 f and I* a r 4 /V 0 . We choose power constant 

r\ 2 

(i.e., f z r Q = constant) since it is quite a safe assumption although in 

practice one would probably shade this with some increase in power. Thus 

2 

if we let V Q = 1/2 Mra u v v /e where v v is the vehicle velocity, and maintain 
Am/m = 1/Mm, T / B = .05 and = 1/2 (1 + T/B) , then 


r 2 M 2 

t+ / \ o m 

I (max) a — — 


2 M 2 

r M 
o m 


(r 


a 


Mv 


5 ^ 

i T - 


+ .00378 


5 M 


+ .00378 


m 


) • 


(90) 


Figure 10 shows plotted as a function of r Q with V Q as a parameter. 

The interesting velocities are 2 km/sec for rockets, 4 km/ sec for "an 
entry mass spectrometer" and 8 km/ sec for earth satellites. 


ij (max) = 8.43 x 10 13 amps at P g = 10 4 torr and T g = 400°K when 
= 50, t = 5 x 10~V 6 - 10 _5 , r 0 = .00254, f = 4 x 10 6 Hz, and V q = 4.3 V 
(M = 50 at 4 km/ sec) . 


It is noted that relatively large instruments are required if the 
sensitivity is to be maintained. Alternative scaling factors would have 
required large powers to accomplish the same end. 


CONCLUSION 


It has been shown that an approximate optimization scheme can be 
applied to a combined quadrupole mass analyzer and a generalized ion source 
to produce the framework for planning specific missions. While numerous 
assumptions and extrapolations have been used, the end results have been 
compared to the experimental data of several instruments and found to be in 
agreement within about a factor of two over a variety of operating conditions. 

Additional effort is required to convert these results to a complete 
instrument for a defined purpose. Size, shape, weight and power considera- 
tions, as well as the nature of sampling and detection of ion current are 
all essential to successful performance. A variety of instrument charac- 
teristics can be shown when these additional factors are related to the 
optimized analyzer. 
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An interesting result has been obtained in that a double-focusing 
magnetic mass spectrometer has a very similar equation to the quadrupole 
analyzer having the form 

k (f ) 3/2 (*? (-i) -%) 

where 3 = Av/v, B = magnetic field r m = magnetic radius and = 
aberrations in a and £5, 
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APPENDIX 1 


QUADRUPOLE MASS SPECTROMETER 
FOR ENERGETIC PARTICLE MEASUREMENTS 

L. G. Hall and M. R. Ruecker 


SPACE SCIENCES LABORATORY 


PERKIN-ELMER 

Aerospace Division 

2855 Metropo litan Plac e 
Pomona, California 91767 



ABSTRACT 


Analytical expressions are developed which define the initial conditions for 
100% transmission of a Quadrupole analyzer and which enable combination with 
equations for the focal properties of an ion source. The equation for the 
combined properties can be used to calculate the ion current intensity of- the 
system as a function of the gas density, parameters of the mass spectrometer 
and upper limits of acceptance of initial energy and initial angle. It is shown 
that the measurement of density is, in principle, independent of the initial 
energy and initial angle at the time of ionization within the established upper 
limits. The equations are particularly useful in design of instruments for 
measurement in the upper atmosphere. 

This work supported by the Jet Propulsion Laboratory and by Goddard Space Flight 
Center of the National Aeronautics and Space Administration. 
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QUADRUPOLE MASS SPECTROMETER 
FOR ENERGETIC PARTICLE MEASUREMENTS 


INTRODUCTION 

In terms of optimum performance, the design of any mass spectrometer considers 
the initial conditions under which the ions enter the analyzer. The initial con- 
dition of the ions, in turn, relates to the ion source properties. It is the 

1 

purpose of this paper to present a means of analytically relating the quadrupole 
analyzer with an ion source in one of two modes of focusing. For this purpose it 
is necessary to provide analytical expressions for both the analyzer and ion 
source operation. 


QUADRUPOLE ANALYZER 

The discussion of the quadrupole analyzer is limited to operation of the quadru- 
pole in conditions of 100% transmission. This choice dictated by a desire to: 
1) obtain a more accurate absolute value for the ion current as related to gas 
density in the ion source; and 2) obtain "flat- topped" ion intensity peaks which 
can be measured by stepping from peak-to-peak. This method of measurement may be 
useful in space applications , where time or information capacity is limited and 
where the peak intensity must be independent of small changes in voltage or fre- 
quency. The discussion that follows thus treats the top width of the peak and 
the top to base ratio as important parameters of the operation. 


The functional Characteristics of the quadrupole analyzer have been adequately 

1 * 2 

described by Paul, Reinhard, V. Zahn , and Brubaker * Their work is not exten- 
sively reviewed here. 


It has been shown that the maximum amplitudes of motion are inversely propor- 
tional to the stability parameters (1 - j8 ) and /3 . Paul, Reinhard, and V Zahn 

x y 

have also computed the stability requirements as a function of entering phase 
angle under conditions 1) that the initial transverse velocity is zero or 2) that 
the initial position is zero. Since here we require an analytical expression 
with all initial conditions (i.e. , position, x q , y , transverse velocity, x q , y Q , 


* References I through 5 herein. See final page of this Appendix. 
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and phase angle , <p) as variables, by retaining only the first few and important 
terms, the following approximate expressions for x and y motion are derived. 


sl "(%- + * x )- (=_ 2 ) x ■ sin [( l - »*]} 


and 


where 


f 


y ) [ c c +2 ( c 2) 


y » A it s±n(-^ + | c_ + 2 (c, ) cos Lt +A 

y 


V " f (* 0 * V*) 

V 6 = f (y 0 , y c , *) 


<t> = Initial phase angle 
C + c 

°0 " 1 ‘ (^) - 2 2 ~ 2 - °- 168 * £ (« y ) 


A1 (1) 


A1 (2) 


°- 8 < N . 


< 1.0 


It is noted in Equation (1) that when p is nearly 1 the arguments — 5 — and 

0 X 1 

1 - - ' 2 " ~ are nearly the same and since the coefficients are of similar magnitude 

one obtains an amplitude modulated wave without carrier which is centered about 

the angular frequency . On the other hand, when f3 is small the y motion is 

^ y 

0yU>t 

the low frequency wave sin — ^ — , with superimposed fine motion (1 + 0.336 cos cot) 
Brubaker has shown the trajectories for both motions computed with specific con- 
ditions and has also derived the elements of Equation (2) by the integration of 
incremental forces on the ion. We differ from his value for 2 A ^ 2 ] only in ob- 

V )y 

taining 0.336 instead of 0.353. 
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From Equations <1) and (2) one can obtain values from the maximum amplitude of 
motion in each plane as a function of initial conditions. Since the y motion is 
simpler we use this for an example as follows: 


From (2 ) 9 


y = A sin 


/'0y“ t \r ,1 

I— ' Xj- + J j 1 + 0.336 cos (ot +<£) 


and we first observe that (y ) = 1.336 A 

•^max y 


Evaluation of (3) and its time derivative at t = 0 gives: 


y = A f sin & 

° y 1 y 


and 


= A 


( f l COS *y " f 2 ** Sin 5 y) 


where 


= (1 + 0.336 cos <f>) 

f 2 = 0.336 sin . 


and 


Solving for A^: 


2 
A 

y 


' ^l y ° 2+ ^ +y ^)( . 


Substituting in (4) gives: 



A1 (3) 
A1 (4) 


A1 (5) 


A1 (6) 


A1 (7) 
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It can be easily shown that for resolution better than = 10 


The term y^ can be neglected. Consequently, 


j | 2.672 

y ° + f2 

PH P y f i 

~ + y o~ 


A1 (8) 


where the absolute signs are obtained in the process of taking the square root. 


By the method employed by Paul an approximate equation for 

lx evaluated at & = 0.92 is: 

max H x 


= 16 


= 0.92 


.35 ^ 


„ XX X 

f x + f. — 7~+ — “ 

3 o 4 <o 5 co 


o 2a 


’ A1 (9) 


where 


= 0.40 - 0.20 cos <p - 0.20 cos 2<p, 


= 0.93 sin <p + 0.26 sin 2 <p uf and 
f,_ = 1.63 + 1.86 cos^+ 0.26 cos 2 cp , 


The variable <u may be eliminated with the following substitutions: 

/ 2ev i \ % . 

y —J 


A1 (10) 


where eV^. is the injection energy and a the entrance half-angle. Using the qua- 
drupole equation: 


2 2 
a co r 
o 


A1 (11) 
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and letting 


y = 



sin a s 


equation (10) becomes: 



In the region of the apex of the stability diagram q = 0.706* 


A1 (12) 


A1 (13) 


Thus dividing y by y in Equation (8) and substituting Equation (13) : 
max o 


^max 2.672 

0.595 r Q y 

y 8 f, 

o M y 1 

y 1 fi 

J o 1 


A1 (14) 


Similar substitutions may be made for x 

max 


From these equations we can now show examples of combined initial position and 

angle in Figure 12 where (x ) and (y ) are plotted as a function of initial 

max max ^ 

phase angle, <^>, with y as a parameter. The value — - 30 is used where 

r g = x q = y Q is the radius of a round entrance aperture; 


Comparison with the curves which Paul has computed for x.y = 0 or x , v =0 

o^o o’ J o 

indicate that the approximate solutions are in good agreement in the x - z plane 
where maxima are within 7% and the same shaped curve is obtained. In the y - z 
plane, ion motion reaches a maxima for a phase somewhat shifted from Paul's, 
however, maxima agree to better than 2%, It appears that one may use the approx- 
imate expressions with less error than 107. in the maxima of the curves. 
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max 


> X max 


By comparing the | and curves it is seen that for equal values of y, 


It should be noted that the maximum amplitude in y - z plane occurs in the useful 
region between 127° > (f 3 > 142° . A study was made over a variety of conditions 
where it was found that for equal initial conditions in the two planes, nominal 
worst case condition is at <f>- — in the y - z plane. At most, a small per- 
centage error will occur, which is in a direction to reduce previous differences. 


Consequently, the following generalizations are made: 


For 100% transmission |y | r 

“I*- 3^/4 “ °» 


thus 


r o _ r o A 2.672 


B f i 
h y 1 


X o y f 2 

0.595 — 2 — + — —■ 
r f, 

e 1 


j 3* 


where 


Thus letting 


w 


- 0.763 


3jt/4 

r 




= 0.238 . 


3rr/4 


o 3 

— and evaluating (14) at $ = gives 

e 4 


_o > 1,09 


fi. 


1 + 1.9 


or solving explicitly for 


in the limiting case: 


fo = 1,09 

r (3 - 2.07 y 
e max 


A1 (15) 


This equation asks: What is the largest value of y entering at any phase and 

for given j3 y which will just reach r Q , the position of the rod. 
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Referring to Figure 13, the interpretation is as follows: As the scan line passes 

through the stability diagram, jQ is at first very small. 


Consequently only low angle ions contained in low values of 



sin a 


r 

are transmitted within a given — ratio, As 

e 

accepted and transmitted by the analyzer with consequent increase in collected 

current until one of two conditions become predominant. If, by some means, Y 

is limited, a region of 100% transmission occurs with a resultant constant 

collected current or 11 flat- topped" peak. If Y is unlimited the x stability takes 

over at a position where j3 = 1-/3 and by the same mechanism ions are lost to the 

y ^ 

x poles as /3 — ►!'. In this condition, a pointed peak is produced. 


increases , larger angles are 


While one could easily plot families of curves relating to specific values of 
r 

— and y and match these to an ion source, this relationship can be obtained 
e 

more explicitly. Therefore, we proceed to examine ion sources. 


ION SOURCES 

When reduced to their elements, most ion sources can be typified as focusing in 
one of two modes. Figure 14 shows the most common method wherein the source ion 
current is peaked up by producing the smallest line or point "image" and thus 
the highest density. However, it is noted that given some thermal distribution 
in the gas being ionized, those with higher side velocity are not transmitted. 

o 

If only a specific side velocity ±y^ is present such as encountered in satellite 
motion the line spectra are not transmitted at all. 

An alternative mode of focusing is shown in Figure 14b where an image of the ion 
object is produced. In this case each position in the image correlates with a 
position in the object. Also, up to a limit determined by the lens aperture 
sizes, the position in the image is independent of initial energies of the ions. 
It is noted then that the distribution of initial energies lies in the angular 
distribution. 
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Since this mode of focusing is more appropriate for energetic particle measure- 
ments, we may proceed to describe the ion density more analytically. For a 


J + 

O 

magnification, M the ion density at the image plane is — r for J current 

FT ° 

per unit area produced is the ionizing beam. Thus the source current is 



A1 (16) 


where r^ is the radius of a round aperture. 

From Abbe's law , the bundle of rays produced within an included angle at a 

point in the object plane, fall within an included angle at the image plane 
by the relationship; 


% 

y V sin 0 
J o o a 


y . V. sin $. , 
3 1 x 1 


A1 (17) 


where eV Q is the initial ion energy and e\L the total ion energy at the image 
plane. 


Letting Q = 


sin 61 


jr_ 

2 


to include a worst case angle and since 




A1 (18) 


From (18), for a given value of (V ) , we note that sin Q. can be kept within 

some upper limit acceptable for transmission by the analyzer by adjusting the 
magnification and thus from (16) the ion source current becomes established. If 
the value chosen for r g is also sufficiently small for 1007. transmission then 
the collected current for the mass spectrometer becomes. 



A1 (19) 
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Consequently for this mode of focusing we find that the measurement of partial 

density is, in principle, independent of initial conditions up to a maximum 

initial energy (eV ) .It remains' to establish r and sin ft* such that 
bJ o max e 1 

the maximum ion current is collected under the restraints of 100% transmission 
and limited initial energy. 

Combining {(19) and (20) yields 

sin 2 ctj A1 (20) 

max 

in which the injection voltage is used for the image plane energy and 

the maximum injection angle o ] . is substituted for 6^. In Equation (20), M 
is an implicit variable used to limit sin 2 a as (V ) is chosen. Since this 
equation contains r and sin a it may be solved simultaneously with the qua- 
drupole equation. 



Combined Quadrupole and Ion Source 

2 

By substituting for y in (15) and solving for V sin a , one obtains 


t- sin 2 a = 0.233 V [/ 

I max ac 


S - 1.09 


? 


A1 (21) 


where a is the largest acceptable angle for given values of B and — 
max y r 


Angles less than this will be transmitted. Further, one may conclude for given 
r 
r 


e ' + 

a and — that the collectable current, I , reaches a peak when a = a x 
I r c max l 


and remain constant when a > a T 

max I 


Substituting in (20) , 


+ 

I 

c 


0.233 V nr 2 
ac e 





A1 (22) 
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+ * e 

where we note that I = o at r = o and at 1.09 — = B . 

C A T r'v 


r <“y 

o J 


These conditions represent zero aperture and zero acceptable angle respectively. 
Zero angle in turn implies infinite magnification and thus zero current density. 

r 

By differentiating for a maximum, one finds that and — contribute equally 

o 

to the "beam width" when 


«) 


- 0.211 ft 


A1 (23) 


and from (21) 


.. . 2 
Vj sin 


0.0583 V B 

ac • 


A1 (24) 


With reference to Figurel3 > for "flat topped" peaks the ion current arrives at a 

■maximum before the scan line reaches the P y - 1 - 0 X intercept. If the 

base width, B, of the peak is defined to be that portion of the scan line within 

the stability diagram (i„e., the theoretical base width); and the top width, T, 

that portion lying between lines of ^3 = 1 - ft (not at the intercept) and 

y x 

where the ion current is at its maximum value, then by calculation of the 


intercepts , 


M - ( 


- - °- 164 ir 


(■ ■ ?)■ 


A1 (25) 


Equation (25) describes the j8 and jS intercepts with the scan line which are 

y dm ^ XT 

required to produce the resolution — for a top ^to base ratio, — . Iso - ~ 
lines may be constructed by inspection. 

We now require that the ion current reach maximum value at the above - scan 


line intercept. This consequently defines the values of 


and sin a for 
max 
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maximum ion current as: 


(Is) - 0.0346 (l - 2) 

max 

and 

V T sin 2 a = 0.00955 V — (l - 

I max, ac m ^ By 

From ''22) , (23) and (26), 


A1 (26) 


A1 (27) 


I + - 3.15 xio- 4 -r, 2 J 0 + (^s) 2 (l-l) 2 . “ C23) 

max o max \ / \ / 

Throughout the T region then I + is constant at this value since 

max 

a max ^ a I * 

Thus, with given limits on the initial conditions and choice of the quadrupole 
parameters one may calculate the expected current of mass peaks meeting criteria 
of shape and resolution. The interesting result is that, for this mode of ion 
source focusing, the collected current is, xn principle, a measure of partial 
density in the ionizing region independent of the initial energies up to . 

An instrument for atmospheric studies has been constructed along these lines 
having r^ = 7.5 mm, r g = 0.25 mm and length £= 22 cm. While the ion 
current exceeds that expected, an accurate confirmation of Equation (28) has not 
been made. This is partially due to the difficulty of obtaining a rectangular 
energy distribution and partially as discussed later to a method of biasing the 
quadrupole rods. 

A check on Equation (27) was made however, to determine the validity of injection 

T 

angle versus resolution and — ratio. Figure 15 shows two partial spectra typical 
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of those used to calculate the apparent injection angle, given measurements for 

T 

the other variables. V ac was the independent variable while — was the pri- 
mary dependent variable. The theoretical base width was used (i.e., without 
tails) measured by extending the peak sides. Figure 16 gives the apparent angle 
vs V where the inset shows the directly measured angular distribution. The 

3.0 

apparent angle is approximately constant over the range of V however, with 

3 C 

the density distribution less than the measured angle by a factor of ~ 1.5. 
This contraction in apparent angle is believed to be due to the effect of qua- 
drupole biasing. 


Quadrupole Biasing 

The length of the quadrupole rods is required to have a certain minimum value 
which depends upon the values of several parameters, von Zahn^ stated this 
relationship in the simple form: 


n 


A 



A1 (29) 


where n = number of cycles of the r.f. voltage which occur while the ion is 
in the quadrupole field. This equation implies that if a half height resolution 
of is required then n cycles must occur in order that the peak tails do 
not make an appreciable contribution to the width of the peak. By relating n 
to the axial ion energy, and using the quadrupole equation one obtains: 


JL. K 

r 

o 


f ac 


Am 


m 


A1 (30) 


where £> is the length of the quadrupole rods, and K is a constant which de- 
pends upon the exact criterion used to measure the resolution and the contribution 
of the tails to the peak base width. If von Zahn's criterion is used, K = 3.5. 

For one percent peak height resolution K « 7. 


In normal quadrupole operation the ion source injection voltage is equal to the 

axial kinetic energy of the ions in the quadrupole, i.e., V = V , When an ion 

X z 

source is required to capture and focus energetic particles which have an 
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appreciable angle of incidence, the injection voltage, V^., must be large. This 
implies that the length of the quadrupole rods will be excessive. It was found, 
however, that by placing a bias voltage Vg between the ion source and the qua- 
drupole rods the mass spectrometer may be successfully operated in a mode where: 

V 2 - Vj - V B A1 (31) 

A schematic representation of this technique, along with an approximate axial 
potential distribution, is shown in Figure 17 >■ The effect of the quadrupole bias 
voltage is shown in Figure 18 . As the bias voltage is raised the peak tails are 
diminished. There is no observable decrease in sensitivity as V„ is increased 
until the point at which the injection energy of the least energetic ions is ex- 
ceeded in which case they are cut off. Even though the reduction in the size 

T 

of the tails is not considered, the — ratio apparently improves as the bias is 
applied. This may indicate that the field around the entrance nozzle is causing 
a contraction of the effective angle of injection. 


APPENDIX 1 REFERENCES 

1. W. Paul, H. P. Reinhard and V. von Zahn; Z. Physik 152 , 143 (1958). 

2. W. M. Brubaker; Neuvieme Colloq. Spectroscopicum Intern., Lyon, 

France, 1961. 

3. Reference 1, Fig, 4. 

4. c.f.: Klemperer; Electron Optics, Cambridge University Press, 1953, p 14. 

5. U. v. Zahn: Diplomarbeit Bonn, 1956 (unpublished). 


59 



THIS PAGE INTENTIONALLY LEFT BLANK 


60 



MHWM& PAGE BLANK NOT FILMED 


APPENDIX 2 
PEAK TAIL EQUATIONS 


In order to obtain a relevant mass peak-width for the quadrupole it it neces- 
sary to describe the "tg.il" distribution as a function of a Aq from the 
center or edge of the theoretical peak. For the 100% transmission portion 
of the peak, it was sufficient to describe the extreme particle as defining 
the top edge of the peak. Since the tail has been shown to be a function 
that continues to infinity, it is necessary to define the tail function be- 
fore a selection of a reasonable cut-off point can be made. 

The purpose here is to obtain a relatively simple analytical expression for 
the tail as a function of Aq t away from the peak in order to handle this 
function in the optimization procedure. Computed data is used to elucidate 
a function that is then matched to experimentally measured peak tails. 

Rather than following a rigorous and ambiguous mathematical derivation, it 
is further the purpose to gain some physical insight into the unstable par- 
ticle behavior. ^ 

We start by replotting a limited amount of computed data from the Mask Study. 
Figures 19 and 20 show the log of normalized particle amplitude, A = y/r Q as 
a function of number of cycles of travel in the quadrupole. Only the enve- 
lope of motion is shown. Phase of particle entry, <J> 0 , and fractional 
distance, Aq^/Aq from the theoretical peak edge are used as parameters as 
well as riy and yy* The data available allow the plotting of y-axis motion. 

An interesting result is found in that after sufficient travel, motion ap- 
proximates exponential growth versus the number of cycles, n. It is seen 
that near entry the motion envelope is non-exponential. It is also observed 
that as Aq t /Aq approaches zero (near the theoretical boundary) the motion 
falls away from the positive exponential growth. This, of course, is to be 
expected since a continuous function must describe the transition from un- 
stable growth to stable, bounded motion. The end function used here will 
then be limited in its applicability to Aqt/ Aq = 0.1, i.e., 10% outside of 
the theoretical boundary. 

We are basically interested in the number of cycles taken for the particle 
to reach the characteristic radius, r 0 , as a function of the initial con- 
ditions. For the purposes here the details of motion at small radius or 
near entry are of no interest if 'we can properly describe the motion as the 
particle assumes a final function. 
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Consequently, we note from Figure 19 and 20 that: 

a. The slope of the positive exponential depends upon 
Aq t /Aq. 

b. The intercept with n = o depends upon <{> 0 and q and y. 

c. Amplitude is proportional to n and y when taken separately (see 
Mask Study) . 

Placing the observations of the steady-state function in mathematical 
notation 


in A = an + f (q,y, <j> Q ) + f(Aq fc /Aq) 

t 

We can either leave f (q,y,<j\j) *- n present position or combine it with A. 
Since we wish to separate the variables and n and y control A linearly, we 
choose to place it in the form 


in * an + f(Aq t /4<l) 

Fitted empirically to the group of curves of q = 

A = f(n,y, <^ o ) exp 

Thus Aq t /Aq not only controls the slope versus n 
the interesting crossover. 

Since the function of n, Y, and <{> 0 appears to be 
particles , we let 

1 t 

f0l>Y 5 <l> o ) = (^Yy + F 2^» 


^ . 043 + .121 (Aq t /Aq)j n 


.0025, we find, 


-1.1 (Aq t /Aq) 


(A2-1) 


but also the intercept and 


of similar form to stable 


(A2-2) 


or 


A = ^ F l Y y + F 2 T1 y^ 6X13 


^.043 + .121 (Aq t /Aq)jn -1.1 Aq t /Aq 


(A2-3) 


i t 

where Fj_ and F£ are assumed to be functions of (f> 0 alone. We now have an 
approximate expression for normalized amplitude that provides linearity as 
a function of q and y individually, and contains a positive exponential 
slope proportional to n. When A = 1 and particles strike the rods, the 
"position", n, is controlled by the initial conditions of q , y and 6^. 
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Since we are interested only in the position at which the particle reaches 
r 0 , we let A = 1 so that the relationship between position and initial con- 
dition becomes , 


( F lY y + F 2 r| y^ ) = ex P 


- ^ .043 + .121 (Aq t /Aq)jn +1.1 Aq^/Aq 


(A2-4) 


Turning now to particle density, for each n> Y and <j^ particle injected into 
the quadrupole there will be a position n of particle loss (either to the 
rods or transmitted to the collector) for each Aq t /Aq assumed. Since we do 
not know the functions anc ^ F 2 we must assume a variety of independent 
values. We further assume they may he negative but have bounded negative 
and positive values. We now let n and y have uniform density distribution 
to +(r|y) and +(yy) where the brackets ( ) indicate maximum values. 

Thus : 


3N 

3tl 


+ (y 

= constant and 


3N 

9y 



+( V 

= constant 



If we integrate the distribution over all n and y values for a given <f> Q we 
can show in general the form to be as follows; 



where we let ft be 


an amplitude function such that, 


ft 



+ F o n 


V* 


3 


(A2-5) 


I ! 

when Fi and IV have specific constant values at cj) 0 = (j) 0 as n and y vary 
» d a 

to change ft. 
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As initial phase changes, now, we expect and F 2 to assume a variety of 
values which are small. Consequently, if we assume 


SM 

n 


2n 

= constant; 
o 


when the value of £2 is small, the constant injection rate creates a high 


particle density within the small band of £2 values. Likewise with large 
Fx, F^, and thus £2, the density is low. These various conditions are 






We further assume that, as in describing the amplitude of stable particles 
(see Appendix 1) , the negative and positive values of ft create envelopes 
of motion on negative and positive sides of the x-axis (in the y plane) . 
Since it makes no difference which side of the axis is involved, we combine 
the negative and positive distribution by letting 


ft 


’ i + F 2V 


(A2-6) 


We also now use absolute values in (A2-3) and (A2-4) . Now let us mentally 
sum all of the combinations envisioned above and assume the following 
distribution: 



This function now represents 


N(ft) = N( Yy ,n y ,<i> o ) 





N(ft) dydrid^ 


(A2-7) 


(A2-8) 


While the result assumed is not a rigorous solution, it provides necessary 
simplicity until a more exact function can be established. It further may 
be a conservative viewpoint since the highest density occurs at low values 
of amplitude where the particles are more freely transmitted by the quad- 
rupole. This will give rise to an assumption of larger tails if the true 
function has a lower relative distribution at low ft. 
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Combining (A2-6) and (A2-4) 


- exp ^-^.043 + .121 (Aq^/Aq) j n + 1.1 Aq t /Aq j 
Also in (A2-7) we let N = o at (1 = ^max or so t ^ iat k = ^o/^max« 


CA2-9) 


N(fi) ® N ( 1 - 


Combining (A2-8) and (A2-9) we have an alternate equation, 

N(n) = N q 1 - ^ - - exp - f .043 + .121 (Aq^/Aq) j n + 1.1 Aq fc /Aq 
max ' ' 


useful between n = £ (ft ) and n - °°. 

max 


We can obtain a picture of the relative location of particles by superim- 
posing the distribution versus ft on a scale of the exponential of -an as 
follows : 



66 



It is now simple to express the number transmitted per cycle N t as 


N, 


= / N(ft) dft = / N(n) dn 

Jo Joo 


(A2-12) 


and the total number of particles per cycle as 


ft 

/•max 

L 


N(ft) dft 


(A2-13) 


Consequently, 


N. = N ft 
t o 


f 1 - tiM 

\ max / 


(A2-14) 


N_ « 1/2 N ft 
X o max 


(A2-15) 


and the tail fraction transmitted. 


0 = N /N m = 


ft 


t' T ft 


2 - 


ft 


max 


ft 

max 


(A2-16) 


Solving for ft/ft^ we obtain 


fi / Q max * (2 - 0 / 2 ), 0/2 


(A2-17) 


for 0 » 0 /4. Since we are interested in 0 « 10 _Z this assumption appears 
valid. 


The second solution provides a physical solution since 0 + 0 as ft = f(n) + 0 
when n *> oo # 


Eliminating ft from (A2-9) and (A2-17) and solving for Aq^ gives 


Aq » Ag 

H t .121 n - 1.1 


tn 


l- 2 

V 


0 ft 


.043 Aq 


max 


.121 n - 1.1 


(A2-18) 


Since the basis of the values in Figures 19 and 20 for Aq^/Aq was at a 
resolution of 100 where Aq = qAm/m = .00706 we substitute this value. We 
also drop the 1.1 from the denominator for simplification because our 
greatest interest is for n > 35 cycles. We further let 2/ft max = k 0 . Thus 


Aq t = 


.0582 

n 



.00251 


(A2-19) 
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or in generalized form 


Mt - f 


- b 


(A2-20) 


We now treat this as a simplified general form of the tail function. Because 
of the similarity of shape (except for magnitude) between x and y tails we 
assume the same form of equation for both. 

At this point, we would normally integrate and average over a range of n in 
order to obtain the effect of the various z velocities from the ion source. 
However, this leads to a more complex expression and becomes unwleldly in 
optimizing the instrument performance. The net result of this exercise is 
to demonstrate that nearly all of the transmitted particles are strongly 
weighted to come from the highest velocity group. We therefore assume that 
the ions arise from a thin plane in the ionizing region that is located 
very close to the position giving maximum velocity (v z ) . The simplified 
equation (A2-20) then applies. 

We further ignore the mild dependence of fimax on (q) and (y) because of the 
very small leverage exerted on the equation compared to the many orders of 
magnitude in a range that we will consider for 0. This would otherwise 
lead only to a very complex result with q and y in both proportional and 
logarithmic form. 

We thus take the equation of (A2-20) and fit it to the experimental peak 
shape shown in Figure 21. This peak on a logarithmic recording appears to 
demonstrate a classical form of the quadrupole peak shape with the relative 
magnitudes of x and y tail clearly shown. The fitted lines of (A2-20) are 
shown in Figure 21. 

From this procedure we obtain; 


and 


Aq 


tx 


0.0297 

n 


In 


0.0258 

0 

x 


0.00128, 



0.0584 o 0.0667 

-IT- — 
y 


0.00250. 


(A2-21) 


If we let *= €L to define the peak base-width to equal transmission on 
either side we find the combined contributions to base-width to be: 


Aq. = Aq. + Aq t 
tx ty 


0.088 


n 


in 


0.048 

6 


- .00378. 


(A2-22) 


It remains now to combine this contribution to the peak width with the 
stable portions of the equation in order to optimize the transmission. 
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APPENDIX 3 


OPTIMIZATION 

We proceed to optimize the collected current by combining ion source and 
mass analyzer equations. In (60) where 



c 


4 K- K„ V 
1 2 ac 



(v ) 3 (n) Cn ) (Y ) (y ) amps, 

^ a y a y 


and (43) where 


(60) 


^ (1 “ T/B) = 2.55 [F 1 (y ) + F 2 (n )] 2 + 1-3 [F^) + F^)] 2 

m y y 


/.088 D .048 \ 

\ 0 / 


v - .00378, 
z 


(43) 


we temporarily assign the following letter symbols to simplify tracking the 
operations : 

D - kl a -T/B), 

m 


„ (.088 /; *048 \ 

p = Nr in — rj v z - - 00378 ’ 


A = 2.55 [F 1 (y ) + F 2 (n 4 )] , 


b = 1.3 [f 3 <y x ) + f 4 (p x )3‘ 


C = 


4 K. K_ V r 
1 2 ac o 


Thus 


U 

B 


A + B + P, or 

1.3 [F 3 (y x ) + F 4 (ti x )] 2 = U - P - A, 
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and 


Y 


x 




Substituting in (60) , 


*1“ 

I = C 
c 


(v z ) 3 (n x )(n y )(v y ) ^[( U - x ! 3 ~ A ) 


1/2 


‘W 


]• 


(A3-1) 


Treating each variable as independent in an equation of constant coef- 
ficients and differentiating ij for a maximum in n , 




dl* ( it _ p _ a X 1 / 2 

I<r> " ° " 'W H ) - F 4 ( V- 


or 


t , 1 / U - P ~ A ^ 

F 4^V 2 \ 1.3 / 


1/2 


(A3- 2) 


Therefore, at maximum (n ) , 


I t = C <0<Y_> 2 v’v 

c z y y / 


1/2 

/ U - P - A \ l( U ~ P ~ A \ 

\ 1.3 ) X 2 \ 1.3 ) 


1/2 


or 


< - ° ( vp 3 <v<y *7i% < B - p - A >- 


(A3- 3) 


+ 

One could now proceed to expand A and maximize I c as a function of riy and 
Yy. However, the complexity tends to cloud an important conclusion. We 
can visualize the result more easily if we note that the optimization 
process could have started with a substitution of ri x into (60). In that 
case we would have found that 


l/o-P-A \ 1/2 
2\ 1.3 / 


vv 


(A3- 4) 


or thus that 


vv ■ w • 


(A3- 5) 


This first important result discloses that at maximum collected current, 
the contributions of angle (y) and position (q) functions to q-space are 
equal. This is true regardless of the eventual values of F 3 , F 4 , (y) and 
(q) which produce an optimum or whether they are mildly coupled functions 
or not. 
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We likewise could have shown 


F,(Y ) - F~(n ) 
1 y 2 y 


_ _1 / U - P - B \ 
“ 2 \ 2.55 ) 


1/2 


(A3- 6) 


demonstrating that the same is true in the y-plane, 


We now back up to simplify visualizing the result and from (A3-5) and (A3-6) 
find that 


(n x ) = C V md V “ F^ <V 
Substituting in (60) , 

+ 3 2 2 ^3 ^1 

r c - c <v <V <v i; ^ 


(A3- 7) 


(A3-8) 


From (43) and the latter substitutions, as well as (A3-5) and (A3-6) , 


U - P - 2.55 [2 ^(Yy)] 2 + 1.3 [2 F 3 (yj 2 


(A3-9) 


Solving for 


, U - P - 2.55 [2 F,(y )]• 

(v x r - y - 1 — 


1.3 x 4 F, 


(A3-10) 


and inserting in (A3 -8) , 


I + = C (v 
c z 


,3 , .2 h( B - r - 2 - 55 t2 W 1 ") 

} 1.3 x 4 F n F„F, 


2 3 4 


_+ 


Differentiating X for a maximum in (y ) , 
c y 


dl 


d(y. 


y = o = <Y y ) 2 (“2 X 2.55 F x [2 F 1 <y y )]j 

+ 2<y )(u - P - 2.55 [2 ^(Yy)]) 


\2 


(A3- 11) 


(A3- 12) 


U 


- P = 2.55 [2 F ] _(y )] 2 + 2.55 F^y ) [2 F^y ) 3 


= 2 x 2.55 [2 F 1 (y )3 i 
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(A3-13) 


2.55 t2 F^y] 2 - 


or 

2. (A3-14) 

y 2 x 2.55 x 4 Y * 

Inserting this result in (A3- 11) yields. 


I + = C (v ) 3 (U - P) 2 /2 

C Z 1.3 x 2.55 x 4 x 8 


or 

, C(v ) 3 (U - P) 2 

T + = 5 

c 1.3 x 2.55 x 64 F FJ J, 

12 3 4 


(A3-15) 


The second important result can now be noted. (A3-10) could have been a 
solution for (y )^. In which case we would have found that 


1.3 12 F 3 (Y X )] 2 - (U l F) , (*3-16) 

or thus that 

2.55 [2 F 1 (y y )3 2 - 1.3 [2 I^Cy^ } 2 . (A3-17) 

Expanding (A3-17) into the original terms, 

2.55 tF 1 (y y ) + F 2 (n y )] 2 » 1.3 [^(y^) + F^)] 2 (A3-18) 

This is the result which was presumed in Equation (38) and was the sub- 
stance of the solutions in that section. It demonstrates that the collected 
current is maximized when the q-space taken in the x and y planes are equal. 


Again it is noted that this is true regardless of the individual optimum 
values of coefficients and variables -and that future work can be conducted 
to optimize the choices of (n) and (y) without revising these results. 


Summarizing the results, so far, we have 
current , 




1 / U - P 

2 V 2.6 


) 


1/2 


found that at maximum collected 


(A3-19) 
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1/2 


and 

¥VVV'i(it) 

It is also apparent from ff (41), that since G- « 2.26 F. and G 0 
(41) is true, i.e. , 1 12 

<V<V = <g 2 k V • 


Expanding P in (A3-15) , we have nox?, 

,3 


+ _ C(vJ 
I c = 212 F, 


( TT ,-.088 o .048 ^ , „„„„„ \ 2 

l D - ( — lri — r ) v z + - 00378 j 


Differentiating I as a function of v 
c 2 


dl 
c 

dv 


+ 3(v )■ 
2 


(2 + 3) 


.088 

ft 

In 

.048 

6 

'•L_L 

,088 


r \ 

ft 

Xn 

j .088 

in 

.048 N 

\ ft 

e , 


.048 


)v d 


0 z 


v + .00378^) 2 

2 1 


v - 3(U + .00378) 
z 


and 


/ 3 (U + .00378) 

(v z> “ 5 


/.088 a »048\ 
\ ft 0 j 


: (f) : 


Cl^-J (U + .00378) 3 (U+ .00378) 2 


212 F 1 F 2 F 3 F 4 


( . 088 d . 
\ ft 


048 \ 


e / 


Letting 

P* 


.088 n .048 
tt d~ 


i = 


27 x 4 C(P + .00378) 

212 x (5 ) 5 (P 1 ) 3 


- 2.26 F 2 , 
(41) 

(A3-20) 
.00378 \ 


(A3-21) 
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Inserting the values for C, U and P and combining the numerical constants, 


I c (max) g 


0.957 K.K 0 V 

1 2 ac 


F 1 F 2 F 3 F 4 V o 


7 (fa) r / . \~ 

— 5 S L. [ (1-T/B) + .00378 ] 5 

( in ^B) 3 W ) 


amps 

(A3-22) 


3 

The dimensions are the same as (60) since (f&) has the dimensions (meter/ 
sec) 3, (A3-22) relates to Mode 3 operation. 


It is desirable to express the component parts of (A3-22) both for later 
calculations and to visualize the relative importance of the parameters. 
Thus substituting P into (A3-19) and referring to (60) for the general 
framework, 


<V 


(n P 


<v 


( V 


(t ) 

y 


-d) 


3 (fJl) 3 (U + .00378) 3 


(.088)' 


1 

2(U + .00378) 

2 F 

5 (2.6) 

1 

2(U + .00378) 

2 *2 

5 (5.1) 

1 

2(11 + .00378) 

2 F 
Z *3 

5 (2.6) 

1 

' 2(U *f .00378) 

2 F x 

5 (5.1) 


} 

1 
1 


3 -3 

in sec 


1/2 


1/2 


1/2 


where 


U = (1 - T/B) 

k m 
m 


(A3- 2 3) 


The third important result which becomes evident after having broadened the 
optimization equation to include an approximate expression for peak width 
due to tails is that the frequency and length have a more powerful effect 
upon instrument sensitivity than other variables. Increasing either of 
these parameters allows an increase in axial velocity without infringing 
the resolution or peak shape. The velocity increase is brought about by 
increase in the repeller-accelerator voltage which allows a faster rate of 
charge formations in the ion source. 
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It would appear from the optimized equation that the resolution could 
approach infinity (TJ -v 0) and still allow a finite current to be collected 
because of the constant, .00378. In practice, of course, this could not 
be true. In applying (A3— 22) the resolution cannot be extended too greatly 
because of the restriction that the base width include a Aq t of at least 
.000706. These limits occur because of the approximate nature of the 
expression for tails. 


Turning now to Mode 1 and 2 (open ion source) operation, from (61) 
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, 4 K. n V r 

V* <0 ( V <V ( 0 <V amps 
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(61) 


which immediately reduces to 
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212 F 1 F 2 F 3 F 4 


(A3-24) 


by Utilizing the same procedures leading to (A3-20) . 
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Now, 


Maximizing I as a function of v by the same approach. 
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Inserting in (A3-24) 
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In this case) 
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APPENDIX 4 


DEFINITIONS 


Peak: The presence or nature of an output (ion current) signal at a 

given mass position. 

Scan: To sweep through the mass spectrum of output signals versus mass. 

Tail: The pseudonym given the nature or "shape" of the output current 

as it decreases away from the edge of a peak. 

Resolution: The ratio, M/AMq, where M is the atomic mass number at 
which a specified peak occurs, and AMo is the width, expressed in 
atomic mass units, of the peak, measured at a specified fraction of 
the peak height. For example, the resolution could be 100 for width 
measured at 10% of peak height. 

Crosstalk: The fractional contribution of signal from a unit peak to 

the adjacent mass positions, in the specified mass range. For example, 
the crosstalk could be given as 10~3 or better in the mass range m/el 
to 100. 

Top-to-Base Ratio: The ratio of the top width to the base width of the 

specified peak, measured at specified fractions of the peak height. 

For example, the top-to-base ratio may be 0.5 or better measured at 95% 
and 5% of the peak height. 

A Flat-Topped-Peak: A peak falling within a specified fraction of the 

maximum peak height over a specified fractional range of the peak base 
width. For example, the peak could be within 5% of the maximum height 
over a range 0.5 of the base width. 

The optimum performance is best described as a peak that is greater 
than a given fraction of the maximum peak height within a specified 
region and less than a given set of values at all positions outside 
that region, i.e., the signal falls within the open region in the 
following figure. 
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Figure for the Definition of Resolution and Peak Shape 
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FIGURE 1. General Configuration of the System Optimized 
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FIGURE 2a. Portion of the Quadrupole Stability Diagram 
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FIGURE 2b. Enlarged Portion of the Quadrupole Stability Diagram 


a = 0.23699 



FIGURE 3. Ion Current Transmitted Through. Quadrupole 
as a Function of q 
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FIGURE 5. Quadrupole Transmission vs Resolution 
Uniform Ion Distribution over n x =* riy - .08 

Y - y = .064 

x y 
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FIGURE 6. I c (max) and Sensitivity 
vs Size (r Q ) 
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FIGURE 7. I c (max) and Sensitivity vs Size (r^) 
With M , P, V as Parameters 
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FIGURE 9 . Sensitivity vs Maximum Mass 
of a Set of Mass Ranges 
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FIGURE 10. Tail Fractional Transmission at the 
Peak Edge vs Rod Length 
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FIGURE 11. Sensitivity vs r at Various V , v and M 
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FIGURE 12. Amplitude vs 
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FIGURE 14. Ion Source Focusing 
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FIGURE 15. Peak Shape 
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TTGURE 18. Effect of the Quadrupole Bias Voltage on Peak Tails 

1-22 CM, r » 7.5 MM, V 1 « 300 V, and V = 400 V 
o J. _ . sc 
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FIGURE 19. Y Axis Quadrupole Ion Amplitude vs Phase Angle of V 
After Entry Into the Quadrupole Field 
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FIGURE 21 . Peak Shape Without Mask 






